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SYNOPSIS 

AUTOMATED OPT ILIUM DESIGN OE REFRIGERATED WAREHOUSES 
AND AIR-CONDITIONED BUIDDINGS 

A Thesis submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 
DOCTOR OF PHILOSOPHY 

fey 

BHAGWAN DAS GUPTA 
to the 

Department of Mechanical Engineering 
Indian Institute of Technology, Kanpur 
August, 1977 


The whole world is striving for the conservation of 
energy at present . All-round efforts are being made to accom- 
plish any specified task with minimum use of energy. Thus, 
there is a need to design systems which are both financially 
feasible and entail minimum functional energy. Most of the 
thermal system designs are based on empirical relations and 
tables available in handbooks and other references which do not 
represent the physical situation truly. The design of refriger 
ated warehouses and air-conditioned buildings falls in the same 
category. 

In the present context the number of refrigerated 
warehouses in the country is multiplying fast and consequently 
the demand on energy is ever increasing . The structural 
cooling load, which depends on the environment, is quite subs- 
tantial apart from the commodity cooling load in such buildings 
The environmental parameters contributing to the structural 
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cooling load are probabilistic and uncontrollable. An unified 
design approach should be thought of for accounting these para- 
meters and reducing the structural cooling load economically. 
Although some work using graphical method of one dimensional 
minimization has been reported in the past, so far no attempt 
has been made in using mathematical programming techniques and 
accounting for the probabilistic nature of the cooling load. 

The aim of this work is to present a unified automated 
procedure for the optimum design of refrigerated warehouses and 
air-conditioned buildings. The design problems are formulated 
as nonlinear mathematical programming problems and solved by 
using multidimensional optimization techniques. 

The thesis has been organized into seven chapters and 
three appendices. Chapter 1 presents a review of the work done 
in the area of optimal design of refrigerated and air-condition- 
ed buildings. The scope of minimum energy and minimum total 
cost designs is discussed. 

Some of the multidimensional optimization techniques 
are stated in Chapter 2 . The interior penalty function 
approach is suggested for transforming a constrained problem 
into a sequence of unconstrained problems . The Davidon- 
Fletcher-Powell variable metric method coupled with the cubic 
interpolation scheme of one dimensional minimization is recomm- 
ended for solving the sequence of unconstrained optimization 
problems . 
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In Chapter 3? the deterministic design of refrigerated 
warehouses with and without a consideration of internal loads 
is presented. A definition of the design day based on the 
hourly outdoor temperature and solar radiations is suggested. 
Several design problems are solved to illustrate the effective- 
ness of the method proposed. The minimization of total cost 
(including equipment, thermal insulation, running and mainte- 
nance costs) is taken as the criterion function. Two types of 
thermal insulating materials and three economics models ar e 
considered to evaluate the total cost. The heat gain components 
are converted into the corresponding cooling load components by • 
using the transfer function method . 

In Chapter 4? the design of air-conditioned buildings 
is considered according to deterministic design philosophy. The 
problem of design of an office building is solved by taking the 
total cost and a weighted average of the cooling and heating 
loads as objectives. 

Chapter 5 deals with the probabilistic design of 
refrigerated warehouses. The probability characteristics of the 
cooling load are derived by using partial derivative method. 
Three different problems are solved by taking the mean value of 
the cooling load and the mean value plus three times standard 
deviation of the cooling load as the objectives. In one problem 
the internal surface temperature is not allowed to exceed a 
prescribed value by more than a specified probability. 
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The application of extremal distributions to the 
design of thermal systems is considered in Chapter 6. Three 
types of extremal distributions are fitted to the hourly maxi- 
mum temperature and solar radiations. It is found that type 
III distribution fits the data most closely. The same distri- 
bution (type III) is found to fit the yearly max imum temper- 
ature data also best. A methodology using the maximum yearly 
temperature data, is developed for the design of refrigerated 
warehouses. The results obtained by using extremal distribu- 
tions and probabilistic design procedure compared well. 

In Chapter 7, general conclusions and contributions 
of this work are summarized along with the recommendations for 
future research in the field. 

The following conclusions are drawn from the present 
investigation: 

(1) It is economical to use thin brick walls with proper 
insulation thickness to reduce the cooling load while 
thicker bricks may be desired from strength point of view. 
Hence a compromise has to be made between cooling load and 
strength in selecting the brick thickness with proper 
insulation. 

(2) Three different orientations of a given building do not 
show significant variation in the total optimal cost; it 
is of the order of 1% while the change in structural load 
is approximately of the order of 2%. 
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The optimization results of the problem in which envelope 
parameters are taken as design variables indicate that, 
for a given volume , maximum possible height is to be 
selected not only for ease of mechanical handling of the 
stored commodity but also for minimizing the total cost. 

An aspect ratio of approximately 1 .5 is found to be econo- 
mical. The cooling load is reduced by approximately 10% 
for such a warehouse compared to the one having a lower 
height with an aspect ratio of 2.0. Same problem is solved 
with the minimization of total cooling load as the objective 
and it is found that the load can further be reduced by 8%. 
Since the desired inside temperature is higher compared to 
refrigerated warehouses , the structural load will be less 
in case of air-conditioned buildings. In these cases it is 
found that it was not possible to reduce the merit function 
by more than 3% of the initial value . It is interesting to 
note that the contribution of the roof to the cooling load, 
is maximu m in any building and hence maximum insulation is 
also desired for the roof. However if one is interested in 
optimizing the total (cooling and heating) load, then the 
insulation thickness may not be maximum for the roof. 
Although the probabilistic design appr®ach required more 
computational effort, it represents a more rational and 
realistic procedure for the design of refrigerated and 
air-conditioned systems. The problem with probabilistic 
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constraint demonstrates a procedure by which it is possible 
to control the inside surface temperature at ary desired 
value which may be required in specific cases. 

In the design example using extremal distributions, it is 
found that the optimal load value is about 4% higher than 
in the case of probabilistic design. This shows that the 
designs based on extremal distributions are more conser- 
vative. 

Finally the interior penalty function method of optimiza- 
tion worked fairly well with all the problems considered 
and is recommended for the optimum design of thermal 
systems . 

The contributions of the present work can be stated as: 
Development of an automated optimization procedure for the 
design of refrigerated warehouses and air-conditioned 
buildings , 

Application of statistical meteorological data in the 
optimum design of such systems, 

Demonstration of the use of extremal distributions in the 
field of thermal systems design, and 
Development of generalized computer programs for the 
opt imum design of refrigerated warehouses and air- 
conditioned buildings. 



CHAPTER 1 


INTRODUCTION 

The ultimate aim of any engineer engaged in the design 
of thermal systems is to produce a financially feasible and 
economically viable design which achieves a desired performance 
level. In the past, when energy costs were low and its avail- 
ability potential was high, the energy required for the 
functioning of any system or the running cost over its lifetime 
was of secondary importance. But now with limited energy 
resources and high costs, it is necessary to conserve energy or 
use it optimally.. A system should be designed for minimum 
functional energy without affecting its performance. In the 
case of refrigerated warehouses or air-conditioned buildings, 
the energy requirements can be reduced by suitably designing 
them. Generally, systems with high initial costs (involving 
sophisticated design concepts) are presumed to have lower 
running costs with high levels of performance, but such systems 
may not be always economically feasible. 

The need for refrigerated warehouses and air- . 
conditioned buildings is increasingly felt in the country. This 
results in the demand of more and more functional energy. Thus 
there is a necessity for improving the design of such systems 

without affecting the total cost much. Generally a number of 

r ■ " .# 1 ' . ' ' ' ■ ' 

alternatives can be thought of for designing a new system or 
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improving an existing one. In refrigerated warehouses the 
structural cooling load, which is dependent on environment, is 
quite substantial apart from the commodity cooling load. The 
environmental parameters contributing to the structural cooling 
load are probabilistic and uncontrollable. An unified automated 
optimum design procedure for refrigerated warehouses and air- 
conditioned buildings is suggested in the present work. A method 
of considering the probabilistic nature of the input parameters 
in the formulation is also presented. The constraints considered 
in the probabilistic formulation are converted into their equi- 
valent deterministic form by making use of probability princi- 
ples. In this work, the systems considered for design are 
confined to buildings only. The design of equipment needed to 
perform the desired function is not considered. 

1.1 Review of Literature 

The work available in the direction of design of 
refrigeration and air-conditioning systems is summarized in the 
following paragraphs . 

1.1.1 Analysis of cooling and heating systems 

In the analysis of cooling or heating systems the 
designer is interested in computing the heating or cooling load 
exactly. The parameters which influence these loads are the 
ambient air temperature, direct and diffused solar radiations, 
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air velocity, characteristic of the enclosure walls and the 
orientation of the enclosure. 'Three methods have been developed 
for accounting the periodicity of external conditions: 
i) Threlkeld's classical approach [ 1 ,2 ,3 ,4] (ii) Transfer 
function method [ 5] and (iii) Finite difference method [ 6 1. A 
computer programme has been developed by Lokmanhekim and 
Henninger [ 7 ] to evaluate the heating or cooling load as well as 
the energy requirement for the enclosure. The programme, in 
addition, may be used for equipment sizing, system simulation 
and economic analysis to provide the desired inside conditions. 
The development of another computer programme, called Gate E, 
which consists of a set of basic programmes in series, is 
reported in references [8,9]. Buffington [10] has presented 
transmission matrix method for the computation of heat gain 
through exterior walls and roof when the outside and inside 
air temperatures are transient. An elaborate design analysis 
of these types is desirable and has become necessary over the 
years due to the ever increasing shortage of energy, rising 
building construction costs and enormous capital and operating 
costs involved in air-conditioning and refrigerating system 
equipment. 

The finite difference method [6] uses an iterative 
procedure for computing the heat transfer through walls. 
However, this technique is not very efficient for composite 
complex structures. In the present work the Threlkeld's 
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classical approach [ 1,2, 3, 4] has been used to compute the heat 
gain due to walls and roof, although this method gives a little 
conservative values of the loads. 

The hourly dry bulb temperature and direct and diffused 
solar radiations data is collected for two cities in India [11, 
12,13 3 • The mean and standard deviations of the above parameters 
and solar-air temperature are computed for all the 365 days of 
the year. This information is further used to select a design 
day for each month. A definition based on the concept of solar- 
air temperature is used to define a design day. Carrier [ 14 J 
gave a definition which accounts for the solar radiation and 
outdoor temperatures separately. A general computer programme 
using the hourly meteorological data is developed to compute 
the heat gain using the classical approach of Threlkeld [ 1 ],. 

In addition to the load due to transient conduction through the 
walls, the programme also evaluates (a) heat gain through 
fenestrations [15 ] (b) heat gain caused by the diffusion of 
moisture through permeable building materials [ 15] , and (c) 
heat generation within the conditioned space that may be due to 
[ 1 5 j (i) lighting [ 14,15,16 j, (ii) human occupation and 
activities [14,15], (iii) electric motors [ 14,1 5] , (iv) 
appliances [14,15], (v) one or more of the processes associated 
with products such as chilling and cooling to the storage 
temperature [14,15,17], (vi) respiration heat released from the 
stored products [15,18] and (vii) infiltration and’ ventilation 
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heat gain [14,15 ,19]. The calculations are carried out by 
using the formulae specified in the cited references and finally, 
the instantaneous heat gain components are transformed to the 
corresponding cooling load components by using transfer 
function method [ 20] . 

Green and Smith [21] have identified the significant 
role of probability theory in the design of building systems and 
associated different probabilities to different events of 
failure. In the present work, a probabilistic optimum design 
procedure is presented by considering the randomness of the 
input parameters. The probabilistic parameters needed for 
evaluating the cooling or heating load are evaluated by using 
the approximate partial derivative method [ 22]. 

Three types of extremal distributions [ 23,24 ] are fitted 
to the daily maximum values of temperature and global solar 
radiation for two cities in India and the best fitted distri- 
bution is used to generate the additional information needed for 
the design. A problem is formulated and solved to show the 
application of extremal distributions in the optimum design of 
refrigeration and air-conditioning systems. 

1.1 .2 Optimization 

With limited availability of useful energy and 
consequent increases in costs, the designers are made to think 
for the optimal use of energy. Serious steps are being taken 
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in this direction through ASHRAE standard 90-75 on energy 
conservation in new building design [25]. Some work has been 
initiated in this direction even earlier. McClure [ 26 ] and 
Stephen [ 27 ] considered the design of building systems to 
reduce the energy requirements. Bonar [28] has analysed the 
different parameters and factors that affect the economics of 
a refrigeration system for a warehouse with reference to the 
minimization of energy costs . 

Stephen [29] has discussed about the economic 
insulation thickness while Spielvogel [ 30 ] has found that more 
insulation can increase the energy consumption. Tobias [31] 
has also discussed some methods for reducing energy requirements 
in residential buildings. Field [32] has given the design of 
a window for minimum energy requirements . Rudoy and Fernando 
[33] have discussed the effect of building envelope parameters 
on energy requirements. Tseng -Yao [34] has described qualita- 
tively the effect of envelope parameters on energy requirements, 
while Ambrose [ 35 ] has studied the. architectural aspects in 
relation to energy conservation. Joseph [36] and Ross [ 37] 
have discussed the recent trends and energy saving procedure in 
the design of refrigerated warehouses. Henize [38] has compared 
the costs of installation and exploitation of high rack const- 
ruction with those of single storied refrigerated warehouses. 
Wood and Bono so [39] have discussed about energy conversion 
efficiencies in different buildings. Stoecker [40] has 
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discussed the different economics models for investment while 
Coad [41 ] has given a methodology for the life cycle cost 
analysis. The first systematic approach for the optimum design 
of air-conditioned systems was presented by Wilson and Templeman 
[42] using geometric programming technique. Jog [43 ] has 
suggested some structural improvements of cold-storage buildings 

1 .2 Objective and Scope of the Present Work 

It can be observed from the available literature that 
the economics part of the design of refrigerated warehouses and 
air-conditioned buildings has not been systematically studied 
so far. Only Wilson and Templeman [42] have made an attempt in 
this direction using geometric programming method . This 
technique (of geometric programming) can only be employed with 
a special class of optimization problems. The one dimensional 
minimization (by graphical method) to obtain the economical 
insulation thickness has been studied by a number of investi- 
gators. However, mathematical programming techniques have not 
been employed to obtain the solution of a general problem so 
far. Also, the design criterion used in conventional designs 
is very much conservative and is based on the peak load value. 
However, this value may or may not occur even once in a year. 
Hence a bore comprehensive definition of the design day over 
a month is used in this work. A weighted average of the 
monthly loads is taken as the yearly load. 
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one to minimize the total cost of the system including the 
initial, running and maintenance costs as well as insulation 
costs and other to minimize the total energy requirements. 

Two economic models i.e. i) Present-worth method and ii) Annual 
cost method are considered to compute the total cost of the 
system. The annual cost model is finally used as the present 
worth method has the inherent drawback of introducing compli- 
cations when the two prospective investments have different 
•lives. However, Coad [4-1 J has preferred to use life cycle cost 
analysis for the optimum investment policy. 

In all the above designs the optimum points are further 
analysed for sensitivity by changing each component of the 
optimum design vector gradually. This study facilitates the 
designer in knowing the parameters that are more important and 
also in knowing the effect of choosing near optimal designs. 
Further, a knowledge of insensitive parameters can be used in 
reducing the number of design variables in subsequent designs. 

1.3 Organization of the Thesis 

The thesis has been organized into seven chapters and 
three appendices. Chapter 1 presents a review of the work done 
in the area of optimal design of refrigeration and air-condition 
ing systems. Some initiative has been taken by the ASHRAE task 
group in finding ways of reducing the energy requirements in new 
building designs. It is becoming conventional now to go for 
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higher and higher heights of the warehouses which help in the 
mechanization needed for the movement and transportation of the 
commodity. Certain values of overall heat transfer coefficients 
are suggested to reduce the energy requirements; hut how to 
obtain these values economically with the existing materials 
still remains a problem. 

Some of the efficient multidimensional optimization 
techniques are stated in Chapter 2. These techniques are well 
developed and have been used in the fields of structural and 
mechanical design. The potentialities of these optimization 
methods have not been fully realized in the field of thermal 
systems design. In the present work, the interior penalty 
function method [44 ] is used for solving both deterministic and 
probabilistic optimization problems. The Davidon-Fletcher- 
Powell variable metric method [45] coupled with the cubic 
interpolation scheme [ 46 ] of one dimensional minimization is 
used to solve the resulting sequence of unconstrained optimi- 
zation problems. Although other methods like Zountendijk' s 
method [ 47 ] of usable feasible directions and gradient projection 
method [48,49] are also available, the interior penalty 
function technique has been selected for the present application 
due to its generality and reliability. 

The deterministic design of refrigerated warehouses 
is considered in Chapter 3- In the present work, the hourly 

at • ■■■■■' . 

cooling load is computed by using the classical approach of 
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Threlkeld [ 1 ]. A definition "based on the concept of solar -air 
temperature is used for the design day and the total meteorolo- 
gical data is processed to select the design day of each month. 

A weighted average of monthly loads is taken as the average load 
over the year. It is assumed that the plant runs twenty four 
hours a day to maintain the desired inside conditions. The 
objective for minimization is taken as the total cost which 
includes the initial, running and maintenance costs of the plant 
and the insulation cost. Two one dimensional problems, are 
involving only a roof and the other involving only a wall, are 
solved. In the first problem, two types of insulation, namely, 
mineral wool and expanded polystyrene are used and it is found 
that mineral wool is economical compared to expanded polystyrene 
Numerical results for the second problem are obtained by consi- 
dering three different brick thicknesses for the walls and it 
is observed that thinner walls with proper insulation are 
economical. A two variable problem is solved to optimize the 
outside surface absorptance for minimum cooling load. The 
effect of orientation on the optimum design of a warehouse (with 
four walls and roof) is also studied by considering three 
different orientations. A more general problem of the design 
of a warehouse is considered by including the envelope para- 
meters as design variables along with the insulation thicknesses 
of walls and roof. The general warehouse design problem is 
also solved with the minimization of the total cooling load as 
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the objective and the two sets of results are compared with 
each other from financial and energy considerations. 

Chapter 4 deals with the automated optimum design of 
air-conditioned buildings according to deterministic design 
philosophy. Here the inputs to the system are the same as in 
the case of refrigerated warehouses but the requirements are 
more rigorous and the system itself may be subjected to a wide 
variation in the cooling load. The inside temperature and 
humidity should be such that most of the personnel inside the 
building feel comfortable throughout the year. This may need 
sometimes cooling and sometimes heating of the building. A 
merit function which takes care of both cooling and heating 
loads, is constructed and the problem is solved for the minimi- 
zation of total cost. The same problem is solved for the 
minimization of total load and the results of the two problems 
are found to be comparable. 

In the Chapter 5, the probabilistic optimum design of 
refrigerated warehouses is considered. Here the design day is 
defined and selected on the basis of mean plus three times the 
standard deviation of solar-air temperature. The expressions 
for mean and standard deviation of cooling load are derived by 
using an approximate method [22] to avoid complicated analysis 
procedure. Three problems are solved to demonstrate the effec- 
tiveness of the probabilistic design procedure presented. In 
the first problem, the merit function is taken as the mean 
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value of the cooling load and in the second it is chosen as the 
mean value plus three standard deviations of the cooling load. 

In the third problem the mean value of the cooling load is 
minimized subjected to a limitation on the inside maximum sur- 
face temperature. In this case, since the constraint itself is 
probabilistic, it is converted into an equivalent deterministic 
constraint before applying the mathematical programming techni- 
ques. This type of restriction has a good physical significance 
if the stored commodity is very sensitive to temperature and 
heat radiation. In fact since the input parameters are probabi- 
listic, this type of design approach is more rational than the 
deterministic one . 

The application of extremal distributions [ 24 ] in the 
design of thermal systems is considered in the Chapter 6. In 
the past extremal distributions have been used in hydrological 
and meteorological problems like floods, droughts, gusts, etc. 
[23]. In a limited sense, these distributions have also been 
used in mechanical design to describe the strength of metals in 
fatigue. They have not been applied in the design of thermal 
systems so far. Three types of distributions for the maximum 
value of outside daily temperature and global solar radiation 
are tried for two different places in India. The same distri- 
butions are fitted to the yearly maximum temperature. These 
types of distributions can also forecast the maximum or minimum 
values at any particular time and some of these forecasts have 
been found to be quite satisfactory. An optimum design 
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procedure using extreme value distributions is presented and is 
illustrated with the problem of design of a refrigerated ware- 
house . 

In Chapter 7, general conclusions and contributions 
of this work are summarized along with recommendations for future 
research in this field . 

Appendix A presents the method of calculating the heat 
gain by the classical approach of Threlkeld [ l]. Appendix B 
gives the data (in tabular form) used for the computation of 
heat gain through different sources, while Appendix C gives 
the tables of coefficients required for conversion of heat 
gain to corresponding cooling load component. 



CHAPTER 2 


OPTIMIZATION PROCEDURE 

Any problem in the area of thermal system design can 
be formulated and solved according to either deterministic or 
probabilistic design . philosophy. If all the parameters affecting 
the design problem are deterministic , the problem can be solved 
according to the deterministic design philosophy. On the other 
hand, if some of the design parameters are random in nature, the 
problem has to be solved according to the probabilistic design 
philosophy. If limitations on the performance and other external 
constraints on the design can be stated and a basis for choosing 
between alternate acceptable designs can be established, it is 
possible to pose the design problem in the form of a mathematical 
programming problem. 

In the deterministic design philosophy a general 
mathematical programming problem can be stated as follows; 

Eind X 

such that f(X) is minimum 

( 2 . 1 ) 

subjected to the constraints 

g (X) i 0 j = 1,. ..,m 

<J 

In the probabilistic design philosophy, however, a general 
mathematical programming problem can be stated as: 
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Eind X 

which minimizes a multivariable function f(X, Y) 
subjected to the probabilistic constraints (22) 

P[gj( X 5 Y) < 0 ] j = 1 , . . . ,m 

where it is the vector of design variables and Y is the vector 
of other parameters affecting the design problem. Here the 
components of X and Y are assumed to be random variables, f 

4 - - 4 - 

represents the objective function and. P[g.(X, Y) <_0J >_p. 

J 3 

denotes that the j"^ constraint has to be satisfied with a 

probability of greater than or equal to some specified quantity, 

p . , j = 1 . ,m , where 0 <_ p . <_ 1 . 

3 i] 

Since f and g. are functions of the random variables 
X and Y, f and g. will also be probabilistic in nature . In 

J 

actual computations, the problem stated in Eq. (2.2) is converted 
into an equivalent deterministic problem to facilitate its 
solution. Eor this the deterministic objective function is 
taken as either the mean value of f or the mean value plus 
constant number of standard deviations of f . The probabilistic 
constraints are also converted into deterministic form by 

making use of probability principles . 

' ' ' ■ # . ' 

2.1 Choice of the Solution Method 

The three general classes of widely used nonlinear 
programming methods are as follows; 
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i) Gradient projection method of Rosen [48 ] which was 
subsequently modified by Goldfarb [ 49 ] ; Though this 
method works well with linear constraints, its efficiency 
is considerably reduced in the case of nonlinear cons- 
traints. 

ii) Feasible direction method of Zountendijk [47]: This 
method is based on the generation of usable feasible 
directions at constraint boundaries . Although this method 
works in a direct manner in solving the problem, the 
analyses during optimization have to be done accurately as 
they influence the rate of convergence and accuracy 
during optimization. 

iii) Penalty function methods: In these methods, the const- 
rained formulation is transformed into a sequence of 
unconstrained optimizations which can be solved without 
much difficulty. These methods are quite reliable and, 
their sequential nature allows a gradual approach to 
criticality of constraints. These methods allow coarse 
approximations to be used during early stages of optimi- 
zation procedure and finer or more accurate approximations 
during the final stages. In the present work, the 
penalty function method of Piacco and McCormick [44 ] is 
used as it has been found to be quite reliable. 

The basic optimization problem of the form (2.1) can be conver- 
ted into an unconstrained minimization problem by constructing 
a function of the form 
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-4 m 

0(X,r k ) = f(X) + r k I G[g.(X)] - (2.3) 

0 1 

where G- is some function of the constraints. If the minimiza- 
tion of the 0-function is repeated for a sequence of values of 
the response factor r k , the solution may be brought to converge 
to that of the original constrained problem given by Eq . (2.1). 


2.2- Eiacco McCormick Interior Penalty Function Method [44] 

In this method, the objective function is augmented 
with a penalty term consisting of the constraints as shown 
below; 


+ + m . 

0(X , r, ) = f(X) - r. I (2.4) 

3=1 g 3 (X.) 

where 0(X, r k ) is called the penalty function, and r k is 

"tld 

called the response or penalty parameter in k minimization 
step. The minimization of 0-function is performed for a 
decreasing sequence of r^ so that 


"k+1 



(2.5) 


The minimization of 0(X, r k ) requires a feasible 
starting point and, in the present work, it is found by trial 
and error. Since each of the designs generated during the 
optimization process lies inside the acceptable design space, 
the method is classified as interior penalty function 


formulation. 
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2.3 Davidon-Pletcher-Powell Variable Metric Unconstrained 
Minimization Method 

In the penalty function formulation, since a sequence 
of unconstrained minimization has to be performed, the selection 
of an efficient method of unconstrained minimization becomes 
very important. All the methods of unconstrained minimization 
find a sequence of improved approximations to the optimum 
according to the iteration 


X i+1 ~ X i + “ S i 


( 2 . 6 ) 


where X^ + ^ is the design vector corresponding to the minimum of 

, -*■•**• 
the 0-function along the current direction S^, is the 

starting design vector and a is the minimizing step length in 

■V 

the direction S^. There are several methods available for 

■f 

finding the search direction in Eq . (2.6). 

In the present work, the David on-Pletcher-P owe 11 
variable metric method [ 45 ] is used. This method can be consi- 
dered as a quasi-Uewton algorithm, and, is a very powerful 

general procedure for finding a local unconstrained minimum of 

th 

a function of many variables [50]. In this method, the i 
search vector S i in Eq. (2.6) is computed as follows 


S, = - [H, ].V0(X, ) 


(2.7) 


where V0(!L) denotes the gradient of the 0-f unction at and 
[H ± ] is a positive definite symmetric matrix. The matrix [Hi] 
is updated according to the following procedure; 
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+ [BL] (2.8) 

(2.9) 

( 2 . 10 ) 

and Y ± = v i^(X i+1 ) - V0(X.) (2.11) 

The updating of [H. ] preserves the symmetric positive defini- 
teness of [H i+1 ] which ensures the stability of the procedure. 

The positive definiteness of [ ] is influenced only by the 

* .. 
accuracy with which <* is determined. To start with [H 0 J[ is 

taken as identity matrix. 


[ H ±+1 ] = [H, ]+ [Mj 


where 


[M.] 


[ N- J 
L 1 J 


* -4- m 

a S . S . 

l l 
+ f? -* 

S. V. 

1 1 


[Hi] V. V. T [hJ 


T 


v i 


2.4 One Dimensional Minimization Method 


In the present work, the one dimensional minimization 

is accomplished by cubic interpolation scheme [46] . In this 

* 

method, the algorithm used to compute <* is reapplied until the 
cosines of the angle between the vectors S and V0 at the 
minimizing step length a is sufficiently small, i.e. 


.i« 


cosO = 


-»• T 

S; 

1 


■ 1+1 


l ? i) • 


< e 


( 2 . 12 ) 
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where e is a small number whose value is taken as 0.05 in the 
present work. This ensures that is minimum along the 

4 - 

direction S^. However, to reduce the computer time, the number 
of cubic interpolations is limited to three. 

2.5 Additional Considerations and Convergence Criteria 

\ 

4 

i) Starting point X q 

* 4 “ 

For the minimization of 0(X, r^ ) , the starting feasible 

• 4 * 

point X q is found by a process of trial and error. Each 
subsequent stage used the solution of the previous stage 
as a starting point. In some cases the overall procedure 
may be accelerated by employing an extrapolation 
technique [5l] to determine starting points for subsequent 
unconstrained minimization cycle. The starting point so 
obtained must be checked to ensure that they satisfy the 
constraints, because at early stage, it is necessary to 
start the unconstrained minimization of 0(X, r^.) from an 
acceptable design point. 

ii) Initial value of r^ 

If r is large, the function is easy to minimize, but 
the minimum may be far from the desired solution to the 
original constrained minimization problem. On the other 
hand, if r is small the function will be hard to minimize. 

■ In the present work the value of r 1 is chosen such that 
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^.5fd 0 ) 1 0(X Q , r.,) < 2 . 0 f(X o ) 


(2.13) 


iii) Subsequent values of 

The total number of r's to be employed is given as an 
input to the problem and the values of r k +1 found by 
using the ratio 


k+1 

r k 


0.1 


(2.14) 


iv) Restarting the [H] mat r lie 

In the case of highly distorted or eccentric functions, 

it might occur after few iterations that S ± . v £>i is 

positive , indicating that \ Is not a direction of descent. 

When this happens, the most efficacious remedy is to set 

[H ] back to [H ] and proceed as if starting over again. 

L i+1 L o 

v) Termination of minimization for each r fe 

For each r k , the minimization of the 0 function is 

terminated when the predicted percentage difference between 

, ,, o\ vaiups is less than a small 

the current and the optimal. 0 -values is 

quant ity , that is , 


<[Hi] v* 4 T v <h < 

= 2 ^; — = ” 2 &T 


(2.15) 


n , ThP value of e used in the 

where c is a small quantity . 

. ~ -„ 1+ n 1 for r, to 0.0005 for 

present woik varied from about 0.1 x k 

higher r fe . 
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vi) Gradient of 0-func.tion 

The gradient of 0-function can he- expressed as 

m 1 

V0 = Vf + r Y . v g, (2.16) 

oil g* 3 . 

However, the gradient of 0-function is computed by using 
forward difference technique in the present work. 


vii) Humber of cubic interpolations 

A maximum of five cubic interpolations are allowed in 
each one dimensional minimization. Out of these, only the 
final interpolation involves the evaluation of the 
gradient V0. All preliminary interpolations use a pertur- 
bation scheme to determine the dot product S . VjZ> as 

S 1 . Vja . (2.17) 

(c + - «-) 


viii) Kuhn-Tucker conditions 

“Hf 

An optimum solution X is characterized by the property 
that there are not other feasible designs X in the 

_ . which corresponds to a 

lower value of the objective function thanf(X ). 
Mathematically this can be expressed by the Kuhn-Tucker 
conditions 


3 f 
8X i 


+ 


l 

je J 




3 8 , 


ax. 

i 


i = 1 ,2, . . . ,n 


0 


(2.18) 
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5. >. 0 , j e J (2.19) 

J 

where %■ axe called the Kuhn-Tucker multipliers and J is 

J - . 

the set of indices of active constraints. Kuhn-Tucker 
conditions are necessary and sufficient for a global 
minimum only in the case of convex programming problems. 
However, they can be used as the necessary conditions to 
test the minimum of any practical design problem. In 
the present work the Kuhn-Tucker conditions are checked 
at the optimum design point. 

ix) Relative minima 

In order to see whether any relative minima exist in 
the design space, two different starting points may be 
used for the sequence of minimization for any example. 

If the two sequences lead to the same optimum design 
(except for small differences that might occur due to 
numerical instability), it can be assumed that local 
optimum is same as the global optimum. 

x) Reducing the total computational time 

It has been observed that some of the automated optimum . 
design problems take a longer time to. satisfy the 
prescribed convergence criteria even after reaching very 
near to the optimum design point. This happens whenever 
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the function being minimized is highly distorted. In 
such cases, it may not be worthwhile to try to reach the 
exact minimum to obtain about 0.5 or 1% decrease in the 
objective at the cost of 40 to 50% more computation time. 
This type of situation can be identified by inspecting 
the progress of the optimization path at various stages. 


% 



CHAPTER 3 


FORMULATION AMD SOLUTION OF DETERMINISTIC DESIGN PROBLEM 


Cooling or heating systems are designed to maintain the 
desired indoor conditions. The outdoor conditions will be quite 
different from the desired indoor conditions. The outdoor 
temperature, solar radiation, air velocity and humidity are 
beyond the control of the designer while the indoor conditions 
are dictated by the nature of the utility for which the system 
is designed. Generally, a refrigerated warehouse, for example, 
is desired for long term storage of frozen or perishable commo- 
dities at a selected low temperature and high relative humidity. 
In such cases the task of the designer becomes extremely * 
difficult in terms of computing optimal design parameters which 
take full care of the uncontrollable outside conditions and 
achieve the desired inside conditions. 

3.1 Solar -Air Temperature and Design Day 

The first problem of the designer is, to choose the 
outside conditions such that the system so designed does not fail 
and at the same time will not be very conservative. To obtain 
the outdoor design conditions, the hourly dry bulb temperature, 
global and diffused solar radiations and wind velocity are to 
be considered at the specified location. Since there is a large 
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variation in the magnitude of these variables over each day of 
a year, the problem is to decide the values of these variables 
based on which the system is to be designed. The general 
practice is to take extreme values- of these variables that 
occurred over a span of years, which are known as design values. 
But a design based on these values will be extremely conservative 
because such design conditions may not repeat even once in the 
total life of the system. In the Carrier handbook of air- 
conditioning system design [ 1 4 2, a design day is defined as the 
day 

(a) On which the dry and wet bulb temperatures are peaking 
s imul t ane ously , 

(b) When there is a little or no haze in the air to reduce the 
solar heat , and 

(c) On which all the internal loads are normal. 

Design Vsed even on the criterion of design day will be over 
safe. Hence a concept based on a combination of temperature and 
solar radiation called solar air temperature , which will cause 
the same rate of heat flow and temperature distribution through 
the material as in the actual situation, is used in the present 
work. Mathematically the solar-air temperature (t ) can be 

t 0 = t Q + (a I /h Q ) + [ eaT^(e'-1 )/h Q ] (3.1) 

where t Q = outdoor air temperature 

= absolute outdoor air temperature 
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a = absorptivity of the surface for solar radiation 
e = emissivity of the surface 
e ' = emissivity of the atmosphere 
cr = Stephen-Boltzmann constant 
h Q = outside surface heat transfer coefficient 
I = total solar radiation (direct and diffuse) incident 
on the surface . 

Usually the third term on the right hand side of Eq. (3.1) is 
ignored while calculating the solar-air temperature as its 
contribution is very small compared to the other two terms. 

Based on the concept of solar -air temperature, a critical day 
is defined in each month as the day on which the average solar- 
air temperature is maximum during that particular month and all 
internal loads remain normal. Twelve such critical days are 
chosen over a year. Thus it is not a single value of the 
temperature that dictates the design but the temperature, solar 
radiation, air velocity and humidity data over twelve critical 
days form the basis for computing the cooling or heating load. 

The computed values of solar-air temperature on the critical 
days are given in Table 3.1 for New Delhi, India. In reality, 
the maximum outside conditions may affect the load after few 
hours of their occurrence depending on the thermal lag of the 
building. By then, the other conditions like internal loads 
may be such that the total load may not be at the peak. Further, 
the peaks occur only for short durations. Hence a weighted 
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average of the cooling load over the twelve critical days is 
obtained, which will be more realistic and safe. 

3.2 Heat Gain, Cooling Load and Heat Extraction Bate [52] 

The terms heat gain, cooling load and heat extraction 
rate as used in the design of cooling and/or heating systems 
are defined as follows . 

3.2.1 Heat gain 

Instantaneous rate of heat gain is the rate at which 
heat enters into or is generated within a space at a given 
instant of time. It is classified in two ways: (i) according 
to the way it enters the space and (ii) according to whether 
it is sensible or latent gain. The first classification is 
necessary to distinguish heat gains due to different modes of 
heat transfer and to calculate energy transfer using appropriate 
equations. The second classification is important for proper 
selection of cooling equipment. 

3.2.2 Cooling load 

Cooling load is the rate at which heat must be removed 
from the space to maintain the desired room air temperature. It 
is extremely important to note that the summation of all heat 
gains at any time does not necessarily equal to the cooling 
load of the space at that instant. The heat gain due to 
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radiation is partially absorbed by the surfaces and contents of 
the space and is not felt by the room air until sometime later. 

3.2.3 Heat extraction rate 

It is the rate at which heat is removed from the 
conditioned space. This rate will be equal to the cooling load 
only when room air temperature is kept constant which rarely 
happens. In most of the cases the control system, in conjun- 
ction with the intermittent operation of the cooling equipment, 
will cause a swing in room temperature. Therefore, the compu- 
tation of heat extraction rate would give a more realistic value 
of the energy removal by the cooling equipment than just the 
instantaneous value of the cooling load. 

3.3 Computation of Heat Cain 

The heat gain from various components of a system can 
be classified as (i) heat gain through exterior walls, roof and 
floor and (ii) heat gain through different sources inside the 
conditioned space. It is evaluated as follows. 

3.3.1 Heat gain through exterior walls and roof 

Heat transfer through walls and roof is computed by 
the classical approach of Threlkeld [ 1 ] given in Appendix A 
for a homogeneous wall. The analysis is extended to composite 
walls by replacing the composite wall by an equivalent 
homogeneous wall . 
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3.3.2 Analysis for composite wall 

The relations for heat gain are given in Appendix A 
for a homogeneous wall. Mackay and Wright [ 4 J have proposed the 
following conversion equations to reduce a composite wall to an 
equivalent homogeneous wall. The ratio of thickness (l) to 
conductivity (K) of the equivalent wall in terms of its values 
for individual layers is given by 


(£) 


eq 


{ K + + ( ^i + ( ^ 2 


4 - 


(3.2) 


where the subscript eq indicates equivalent wall, and i, o, ml, 
m2,... denote inside, outside, first middle, second middle,... 
layer of the compositive wall respectively. 

Similarly the product (KPC) for the equivalent homo- 
geneous wall is given by 


(KpC) = -f- [1.1(6) (KPC) + 1.1 (|) (KPC) + 

eq (1 } K ± i K m1 mi 


\K‘ 


eq 


1 .1 ( K ) (K p C) 2 + • • • ] + 

m2 


(KPC) T 
LQ r(k) 

(i?) 0 

K eq 


0 . 1 ( 1 ) - 0 . 1 ( 6 ) 

K ml m2 


(3.3) 


where p is the density, 

and C is the specific heat of the material of different layers 


of the wall. 



32 


The products S L and S ^Kcan be evaluated for the 
composite wall using Eq. (A..1l)of Appendix kt 


SI = 
n 




/ J n T 2 \ 
1 > 


1/2 


(3.4) 


K 

Substituting a = ^ in Eq . (3.4), 

to y, 0 \ to j2 l/2 

s n L = (if c01 ) = <-r ■ p ■ K[>0) (3 - 5) 

Equation (3.5) for an equivalent homogeneous wall can be written 
as 

(Sn^eq = (? (3 ' 6) 

u ^ eq 


Substituting the values of (^) and (KPC) fromEqs. (3*2) 

K eq eq 

and (3.3) into Eq. ( 3 . 6 ), the expression for (S n ^) e( ^ can be 
obtained. Similarly since 

(S n K) = (|fK 2 ) l/2 (3.7) 

(S n K) e can be evaluated as 

1 /2 

(V^eq - ( (KpO) eq ) (3.8) 

The values of (6) , (KpC) , (S I) and (S I) are to be 

substituted in place of their corresponding quantities in the 
relations given in Appendix A to evaluate the heat gain through 
the composite wall. 
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3.5*3 Solar heat gain 

The incident solar radiation (I) needed for the 
computation of solar-air temperature in Eq. (A. 4) of Appendix A 
calculated by using the hourly data of global and diffused solar 
radiations [12] according to the following procedure: 

Known data: 

latitude of the place = 1 

cl * 

Direction cosines of normal to the surface = (l,m,n) 
Global radiation = 1^ 

Diffused radiation = 1^ 

The values of declination (6) and equation of time (EQT) 

s t 

needed in this computation are given for 21 day of each month 
in reference [15 ]. The values for the other days of the month 
are linearly interpolated in the present work. The values of 

g-f 

declination and equation of time for 21 day of each month are 
given in Tables B.1 and B.2 of Appendix B. 

The hour angle is given by 

H = 0.25 x MIN (3-9) 

where MID = time in minutes from solar noon. It is obtained by 
adding the value of equation of time EQT corresponding to the 
day D of the month to the local time. Then this time is to be 
evaluated with reference to solar noon in minutes. 
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Direction cosines of the solar beam are given by 


COS 

z 

= sin 1 sin 6 
a 

+ cos L cos <5 cos H 
a 

(3.10) 

cos 

w 

= cos 6 sin H 


(3.11 ) 

cos 

s 

= + (1 - cos 2 W 

- co S 2 Z)^ 2 

(3.12) 


where cos S >0, if cos H > tan<5/tan L . 

cl 

The cosine of the incident angle © is given by 

cos© = 1 cos Z + m cos W + n cos <5 (3.13) 

Total solar radiation (I) incident on the wall per unit area is 
given by 

I = (I T - I D ) cos© + I D (3.14) 

Dor horizontal surfaces the values of 1^ are taken 
same as the values reported by meteorological centres [12] while 
for vertical surfaces they are assumed to be equal to half the 
values given for horizontal surfaces [53]. The value of heat 
gain due to solar radiation (I) so obtained is used further in 
the calculation of heat gain from exterior walls and roof. 

3.3.4 Effect of corners 

All the warehouses and air-conditioned buildings are 
three dimensional structures. The exact amount of heat gain 
in such structures can be computed if the effects of corners, 
depressions in walls and other departures from one dimensional 
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systems are considered. As one dimensional approximation 
neglects the energy conducted through the coruers, any estimate 
of heat transfer based on one dimensional approximation and 
inside area of the wall underestimates the true heat gain at 
any given instant. On the other hand, if buildings are treated 
as one dimensional structures and the external wall surface 
area is used in computing the conducted heat energy then the 
true energy conducted will be overestimated. In fact, a two 
dimensional transient conduction equation should be considered 
to account for the corner effects. Jacob and Hawkins [54] have 
given some shape factors to this effect in steady state heat 
conduction. In the present work the external surface area is 
considered for the computation of heat gain through walls, as 
the solution of two dimensional transient equation is compli- 
cated and the shape factors are available only for steady state. 

3.3.5 Outside and inside film coefficients 

The outside and inside film coefficients needed in the 
computation of heat gain are calculated by using the following 
equations [ 6 ]. 

Outside film coefficient: The outside film coefficient, h , 
depends on the wind velocity and the type of surface. For 
different types of surfaces the equations for h Q are given 
below: 
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(i) For very smooth surfaces like glass 

h Q = 3.41 + 0.85 V (3.15) 


(ii) For smooth surfaces such as planed wood and plaster 

h Q = 4.4 + 0.91 V (3. IB) 

(iii) For moderately rough surfaces like finished concrete 

h Q = 6.34 + 1 .21 V (3.17) 

(iv) For rough surfaces such as stucco 

h Q = 6.84 + 1 .52 V (3.18) 

where V is the wind velocity in km/hr and h Q is the outside 
film coefficient in kcal/hr-m -°C. 

In the present work, the value of h Q is obtained for 
the design day by taking V as the average wind velocity over 
the day and the external surface as smooth. 


Inside film coefficients The inside film coefficient depends 
on convection and radiation. The convective component of film 
coefficient, h , is a function of the temperature drop across 

V • 

x 

the film At and is given by the following equations s 
(i) For walls in winter or summer 



0.25 


(0.274 x 4.88) ( At) 


(3.19) 
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(ii) For floor in winter and ceiling in summer 

h = (0.126 x 4.88) ( At) 0,25 (3.20) 

V * 

i 

(iii) For ceiling in winter and floor in summer 

h = (0.295 x 4.88) (.At) 0 * 25 (3-21) 

i 

p 

where h = inside convective film coefficient in kcal/hr-m -°C 
c i 

and At = temperature drop across the film, °C . 

The inside film coefficient corresponding to radiation 
component, h^ , is given by 

e.e 0 T. + At ^ T ^ 

h r _ = (4.88 x 0.1714) ( —“q— ) ~ (t§§) (3.22) 

where e^ = emissivity of the transmitting surface 
e 2 = emissivity of the surrounding surface 
At = temperature drop across the film in °C and 
T^ g = absolute temperature of inside surface in °K. 

The film coefficients h and h cannot be calculated f rom Sqs . 

c i r i 

(3 .19) to (3.22) as they are functions of the temperature drop 

At across the film, while At itself is dependent on h and h . 

1 x i 

Hence a trial and error procedure has to be adopted using a 

standard value of the combined (radiation plus convection) film 

coefficient and then computing At. This approximate At is then 

used in Eqs. (3.19) to (3.22) to obtain the corresponding 

values of h and h . The later values are then used to 
c i r i 

obtain a second approximation to At, which would give more 
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accurate values of h and h . This process is continued till 

i-' * Jl . 

i x 

the values of h and h in two successive iterations do not 

w - X * 

1 1 

differ by more than 0.5%. Then the total inside film coeffi- 
cient h. will be equal to the sum of h and h . However, the 
l c^ r^ 

inside film coefficient in the present work is assumed to be a 
constant since it is not possible to evaluate the temperature 

drop across the film, fit, using Threlkeld method [1 ] for heat 

0 

gain. 


3.4 Heat Gain Through Sources Within the Conditioned Space 

The various factors contributing to the heat gain 
inside the conditioned space are air changes, people, lights, 
equipment products. These are discussed below. 

3.4.1 Air changes 

Each time the door is open, Some outside air enters 
the storage room. The temperature of this warm air must be 
reduced to the storage room temperature which adds to the 
heat gain. The traffic in a refrigerated room usually varies 
with its size or volume. The number of times the doors are 
opened is dependent upon the volume rather than the number of 
doors. The values of average air changes per twenty four 
hours for storage rooms due to door openings and infiltration 
are reproduced in Table 33.3 of Appendix B from reference [15]: 
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p x V x (q - q. ) 

The infiltration load = P~x~T (3-23) 

i 

where p = pressure of inside air 

Y = volume of the air = number of changes per twenty 

four hours x volume of the ware-house 

q - heat associated with outside air at an absolute 
^■o 

temperature of T q 

q^ = heat associated with inside air at the absolute 
temperature T^ and 
R = universal gas constant. 

3.4.2 People 

The rate at which heat and moisture are given off by 
human beings under different states of activity is given in 
reference [14,15]. In many applications these sensible and 
latent heat gains constitute a large fraction of the total load. 
Appreciable variations in the heat emission rates must be 
recognized according to the age and sex of the individual, state 
of activity, environmental influences and the duration of 
occupation. In a warehouse the people are assumed to do 
moderately heavy work during loading and unloading. The corres- 
ponding value of the heat gain is taken from Table B.4 of 
Appendix B. 
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3.4.3 lights 

The instantaneous rate of heat gain from electrical 
lighting, is computed from the following equation [14, 1 5', 16] 

q g 2 _ = Total light wattage x use factor x special allowance 

factor x ~(^ 2 ' " 2046 " x 1 8) ' kcal / hr (3.24) 

The total light wattage is obtained from the ratings of all 
fixtures installed. The use factor in Eq . (3.24) is the ratio 
of wattage in use, for the conditions under which the load 
estimate is being made, to the installed wattage. The special 
allowance factor is introduced to take care of fluorescent 
fixtures and for fixtures which are either ventilated or 
installed so that only part of their heat goes to the condit- 
ioned space. Eor fluorescent fixtures, the special allowance 
factor is taken as 1 .20 in order to allow for the power consumed 
in the ballast. Eor ventilated fixtures, recessed fixtures 
and the like, manufacturer's data must be sought to establish 
the fraction of the total wattage which is expected to enter 
the cooled space. 

3.4.4 Equipments[ 1 5 ] 

When equipment of any sort is operated within the 
conditioned space by electric motors, the heat equivalent of 
this operation must be considered in the calculation of heat 
gain. The equation for calculating this heat gain is given by 



Horse power rating; x load factor x 2545 
Motor efficiency x (2.2046 x 1.8) 


(3.25) 


L em 


where it is assumed that both the motor and the driver equipment 
are within the conditioned space. If the motor is outside the 
space than the factor "motor efficiency" will not appear in 
Eq. (3.25). The load factor is merely the fraction of the 
rated load which is delivered under the conditions in which the 
cooling load estimate is made. The motor efficiency may he 
approximated as 50 to 60 percent at l/8 horse-power rating, 
increasing to 80 percent at 1 .0 horse -power and 88 percent at 
10 horse-power and above as given in reference [15]. 

3.4.5 Products 

A product placed in a refrigerated warehouse at a 
temperature higher than the storage temperature will lose heat 
until it reaches the storage temperature. The quantity of heat 
removed may be calculated from a knowledge of the product, 
including the state while entering the cold room, final state, 
weight, specific heat above and below freezing, freezing 
temperature and latent heat of the product. When the product 
is cooled from one state and temperature to another state and 
temperature, the following equations can be used for calculating 
the heat gain. 

Heat removed, q i , for cooling from an initial temper- 
ature t.j to some lower temperature t 2 above freezing; 
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q i = m C i (t - t 2 ) 


(3.26) 


Heat removed, q^, from an initial temperature t^ to 
the freezing point, t^ , of the product: 


q^ = m x C x (t^ - t^) 


(3.27) 


Heat removed, q^, while freezing the product: 


q ±f = mh. f 


(3.28) 


Heat removed, q„^, from the freezing point to the 
final temperature t^ below freezing: 


q f f = ix C. z (t f - t^ ) (3.29) 

where m is the mass of the product, C and Ch are the specific 
heats of the product above and below the freezing point, t^ is 
the freezing temperature and tu ^ is the latent heat of freezing. 

3 . 4.6 Heat of respiration 


All fruits and vegetables are alive 5 they continue to 
undergo changes during storage. The most important of these 
changes is produced by respiration. During this -process energy 
is released in the form of heat. The amount of heat liberated 
varies with the type and temperature of the product. The colder 
is the product, the lesser will be the heat of respiration. 

The rate of evolution of heat of the product' is given in Table 
B.5 of Appendix B taken from reference [15]. 
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3.5 Rate of .Loading and Unloading 

It is important to consider the rate of loading and 
unloading the warehouse in computing the heat gain. It is 
assumed that the warehouse will be loaded to its full capacity 
in one month, with twelve working hours per day, as shown in . 
Figure 3.2(a).' /v ’ Then the average mass of the product loaded per 
hour will be 

* 1 

m h = 2 z num1:>e:r hours of loading 

average rate of loading (mass) 
number of hours per day 


1 12 x m m 

“ 2 24 ~ 4 


(3-30) 


where m^ is the average mass of the commodity per hour and m 
is the actual rate of loading (mass per hour). 

It is assumed that loading is done only for twelve hours 
a day. Thus the sensible heat load is given by 


q g = m^ x C x (t 1 - t 2 ) (3.31) 

This sensible heat load will be only for the period of loading. 
It is assumed 'chat in the next twelve hour of the day . the 
temperature of the commodity will come down to the storage 
temperature . 

*This model of loading is assumed as most of the refrigerated 
warehouses in India are loaded once to their full capacity with 
the seasonal commodity (like potato) and then the commodity is 
preserved for a long time. However other models like continuous 
loading and unloading can also be assumed if needed. 
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Respiration heat gain: For the respiration energy, the weight 
of the product m r is computed as 

m r = (|r + mr) + (m 2t + ~ + At) + ... + [4(n-l) + 3] 

= [3+7 + • •• + 4 (n-1 ) + 3 ] 

= I 1 . § [6 + (n-1 )4 J = I 1 . f(4n + 2) (3.32) 


where n is the number of days of loading and t is the number 
of hours of loading per day, taken as equal to 12 in this case. 

Thus the average mass of the product per hour will be 


*r 


■f 12 . §(4n + 
n x 24 



f(2n + 1) 


(3.33) 


The mass per hour m r is assumed to be valid only during loading 
period. Once the warehouse is loaded to its full capacity, 
then the mass per hour will be equal to the total mass of the 
commodity stored. When all the above-mentioned heat gains are 
added, one gets the total heat gain. This heat gain is to be 
converted to the corresponding cooling load as outlined below. 


3.6 Cooling load 

The hourly heat gain is to be converted into the 
corresponding cooling load as the radiations absorbed by the 
interior furnishings and the structure reach the conditioning 
equipment after a considerable delay in time. There are three 
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methods available for this purposes (i) Exact method [55], (ii) 
Approximate method known as weighting factor method [ 15 , 52 ] and 
(iii) Another approximate method [ 5 ] known as the transfer 
function method. 

3.6.1 Exact method 


This method was given by Kududa and Powell [ 55 ] to 
convert the heat gain into cooling load. It requires a laborious 
solution of energy balance equations involving room air, 
surrounding walls, infiltrating and ventilation air and internal 
energy sources . The principle of calculation is based on 
considering a fictitious space that is enclosed by four window- 
less walls, a ceiling and floor and having no infiltration or 
ventilation air and internal energy sources. The six equations 
that govern energy exchange at each of the inside surfaces at 
any time are given by 


6 

q . = h (t . - t ) + T g . . (t . 

1,T ' °i 1>T a ’ T As 13 1)T 


t . ) 

H? T 


for 3 = 1 , ,6 


(3.34) 


where 


q. = rate of heat conducted out of surface i at the 
inside surface at time ? 


h = convective heat transfer coefficient at the 
c - 



interior surface i 

= radiation heat transfer factor between the interior 


surfaces i and 3 
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t a ^ = inside air temperature at time t 
t^ T = uniform temperature of the interior surface i 
at time t 

t . = uniform temperature of the interior surface 3 

3 > t 

at time t 

The equations governing conduction within the six slabs cannot 
be solved independently of the above equations since the energy 
exchanges occurring within the room affect the inside surface 
conditions which in turn affect the internal conditions. 
Consequently one is faced with the problem of solving six 
equations simultaneously along with the governing equations of 


conduction in six slabs in order to calculate the cooling load 
at any time t (Q^ t ). 0 ^ t can be expressed as 



A rigorous approach such as this for calculating the cooling 
load would be practically impossible, especially when it has to 
be incorporated in an automated optimum design programme. 


3.6.2 Weighting factor method 

In order to simplify the calculation procedure for 
obtaining the cooling load, various approximate methods have 
been devised to take the instantaneous heat gain components, 
convert them to load components and add these components 
together to obtain the instantaneous load. Generally these 
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methods assume a certain fraction of the heat gain to be 
convective and/or latent (and consequently felt by the space air 
instantaneously) and the remaining fraction to be radiative 
which would not be felt until sometime later (after it has been 
absorbed by the interior surfaces and subsequently convected to 
the air) . The general name associated with these methods is 
"weighting factor methods"'. The simplest one is given in 
Table C.1 of Appendix C, wherein all the components of heat 
gain that normally occur in a space are listed along with the 
percentage of that gain which can be considered radiant, 
convective and latent. The cooling load component corresponding 
to the heat gain is then given by the following equation 


Total value of heat gain 


component of gain which 


component at the time 

X 

is convective plus latent 

+ 

under consideration 


(expressed as a fraction) _ 



Average value of heat 


% of component gain 

gain component over 3? 5 

X 

which is radiant 

or 7 hour period upto and 

(expressed as a 

including the time of 


fraction) 

..calculation _ 


— — 


(3.36) 


3.6.3 Transfer function method 

Mitalas [5] has given the transfer function method for 
obtaining cooling load from heat gain. This method is nothing 
more than a modified -weighting factor method for obtaining the 
cooling load from heat gain. A completely different set of 
factors (or coefficients of room transfer function) are given 
for different types of structures (like light, medium and 
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heavy structures) and in addition, different weights are given 
to different hourly heat gain values included in the averaging 
scheme. The heat gain q p T is given in the form of a time 
series, which indicates the values of heat gain at equally 
spaced points in time. The corresponding cooling load Q T at any 
time v can be related to the current value of heat gain q 

B f T 

and the preceeding values of cooling load and heat gain as 


Qt = J-, (V o 4 e,x + Y 1«e,T-A + - 


( w 1 Q t-a + w 2 Q t . 24 + w 3 « t . 3 a •••) 


(3.37) 


where i denotes the number of heat gain component and A repre- 
sents the time interval. The quantities v Q , v^ , ..., w^, w^, ... 
indicate the coefficient of transfer function given in Table C.2 
of Appendix C. 

In the present work the hourly heat gain components 
are converted to the corresponding cooling load components 
using the transfer function method. 


3.6.4 Cooling load computation 

The cooling load is calculated for all the twenty four 
hours of the critical day for each month. The load, during an 
hour is taken as the arithmetic mean of the six values descen- 
ding from the maximum value. The actual load is taken as the 
weighted average of the twelve hourly load values calculated, 
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one corresponding to each month of the year. The weights for 
different months are computed as follows: 

If t denotes the maximum average monthly solar -air temperature 
e i 

and t^ the inside temperature, the weight for the month i, W , 
is computed as 

W. = (t -t.) (3.38) 

e i 1 

These weights are normalized as 

W . 

W in = iT for i = 1 ,.. .,12 (3.39) 

where t = maximum average solar-air temperature in the 
e i 

month i 

t^ = inside dry bull temperature 

= normalized weight corresponding to the month i 
12 

W = l W . 

i=1 1 

Thus the total hourly cooling load, Q, valid over the year is 
computed as 

Q = J* W in Q i (3 ' 40) 

where is the hourly cooling load for the month i. 

3.7 Economics Models [40] 

The following four economics models are generally 
considered in the analysis of any investment policy. Any of 



50 


them can be used to indicate the most favourable choice amongst 
several alternatives. Some of their merits and demerits are 
discussed briefly in the following paragraphs. 

3.7.1 Pres ent-worth method 

In this method all the costs and incomes are translated 
into the present worth. In order to determine the present worth 
of a future sum, first the interest rate must be established. 
Then present worth is the value of a sum of money at the present 
time that, with compounded interest, will have a specified value 
at a certain time in future. Thus 


P = - — — = (PTO)S (3.41) 

(1 + r) n 


where P = principal 

r = rate of interest per period 
S = total amount to be repayed at future time 
n = number of periods and 
PWP = single -payment present worth factor. 

The series present worth is given by 


P 


(1 + r) n - 1 
. r(l + r) n 


where P = original amount 


(R) (SPWP ) 


(3.42) 


R = amount withdrawn at the end of first and subsequent 
periods and 


SPWF = series present worth factor . 
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Although the present -worth method is widely used, it presents 
some difficulty if lives of two possible investments are 
different . 

3.7.2 Annual cost method 

It translates all non annual costs to an annual cost 
basis . In this method also the rate of interest must be 
established first . 

0 

R = P ] = P(CRP) (3.43) 

_(l+r) n -1 

where CEP = the annual recoveiy factor. 

The annual cost method has two advantages over the present worth 
method. One is that there is no complication introduced when the 
prospective investments have different lives. The second is 
that it is more natural for most people to think in terms of an 
annual cost than in terms of the present worth. 

3.7.3 Rate of return method 

In this method no assumption is made regarding the 
interest rate . Instead the money for the investment is consi- 
dered to be in hand, and the rate of return is calculated as 
though it were an interest rate received by making an outside 



investment . 
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Here we find r such that 


P(CRE) - Salvage value (SHE) = Net income 


where SEE = sinking fund factor = 



(3.44) 


3.7.4 Break-even point 

This method assumes that money is borrowed at a 
specific rate of interest and that the loan is paid off as 
rapidly as possible with no profits extracted. The break-even 
point is defined as the time where the loan is paid off and 
profits begin following to the investor. In contrast to the 
rate of return method where the interest rate was unknown, the 
life is unknown in this method. Thus one finds n such that 

P(CRE) - Salvage value (SEE) = Net income (3.45) 


3.7.5 Annual, cost with simple interest 

A simple economics model of annual cost with simple 
interest instead of compound interest is also considered in few 
cases. Here 


R = 


P(1 + n:-r) 


n 


(3.46) 


In the present work, the present worth and the annual 
cost with compound interest and simple interest are considered 
to decide about the investment policy for some simple one 
variable and two variable optimization problems. However, only 
the annual cost method with compound interest is used in 
further analysis. 
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3.8 Thermal Insulation [15,56] 

The requirement of using thermal insulation in refri- 
gerated warehouses or air-conditioned buildings is a result of 
simple economics principle, that is, to minimize the cost of 
equipment necessary to provide' the desired cooling load. In 
addition, thermal insulation is required to prevent surface 
condensation. The choice of an insulating material usually 
involves a compromise with regard to several desirable properties 
like low thermal conductivity, low heat capacity and, sometimes, 
strength. Consideration should also be given to the properties 
of insulation related to the health and safety including its 
inability to support vermin and insects, freedom from fire 
hazards, dust or objectionable vapours, loose particles that 
may be irritating to the skin, odour, resistance to decay and 
resistance to odour absorption and retention. Moisture should 
be kept away from insulation, to avoid deleterious effects. 

The factors like ease of installation, although, not a 
physical property, is so important that it cannot be over 
stressed. Generally the designer should select those thermal 
insulating materials which can be properly installed with a 
minimum of skill. The expanded polystyrene and mineral wool 
are generally used for low temperature and long life applications. 
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3.9 Water Vapour Barriers [15] 

There are materials which provide resistance to the 
transmission of water vapour under specified conditions. They 
are classified as structural membrane or coating barriers. A 
water vapour barrier does not necessarily stop the flow of 
vapour, but serves as a medium of control to reduce the rate and 
volume of flow. These water vapour barriers prevent the 
accumulation of water within the insulation or construction. 

By retarding the transmission of water vapour, water vapour 
barriers help to (a) keep the insulation dry, and reduce the 
heat load requirements for the cooling system, (b) prevent 
structural damage by rot, corrosion or expansion effect of 
freezing water, and (c) reduce paint problems. 

The primary source of moisture is outdoor air which 
enters while opening and closing the doors, cracks and joints 
in the outdoor structure and discontinuities in the vapour 
barrier or insulation systems . Many construction materials like 
mortar, wood and brick are also permeable to air, which can pass 
through walls when high velocity which impinges on exterior 
surfaces causing partial vacuums in other portions. In fact 
the outer surface of any structure is never air tight. 

Bor insulated construction in which the vapour flow 
may be in either direction for an extended period of time, the 
insulation must be protected by sealing inside with low 
permeance barrier. The barrier should be provided with no 
openings and lap joints as they reduce its effectiveness. 
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The transmission of water vapour through materials under the 
influence of a water vapour pressure gradient is called "vapour 
diffusion". The design equation commonly used for calculating 
the quantity of water vapour transmitted is given by [15] 

W = uAt AE (3.47) 

where W = total weight of vapour transmitted, 
y = permeability, 

A = area of cross section of the flow path, 
t = time during which the transmission occurred, 

Ap = difference of vapour pressure between the ends of 
t he f low p at h , and 
1 = length of the flow path . 

Whenever it is convenient to deal with material 
thicknesses other than unit thickness to which y refers in Eq. 
(3.47), use is made of permeance coefficient M, where M = ^ . 

The corresponding flow equation is given by 

W = MAt Ap (3.48) 

The resistance to water vapour flow through a material 
is the reciprocal of its permeance and is equal to the thickness 
divided by its permeability. The overall resistance of an 
insulated section to the flow of water vapour is the sum of 
resistances of the various materials in series comprising the 
section. The forces other than the pressure of water vapour, 
such as hydraulic pressure , absorption, adsorption, 
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hygroscopic ity and capillarity, affect the water vapour flow 
and cause moisture migration. These various forces are inter- 
related hut the actual relationship is very complex and has not 
been expressed in a simple way for use in designs . 

In many cases, the consideration of water vapour 
pressure difference alone can be used to obtain a fairly accurate 
calculation of water vapour flow. This consideration is useful 
in designing vapour barrier systems or in preventing the possi- 
bility of condensation within a given insulated system. 


3.10 Problem Formulation 


The problem of optimum design of refrigerated ware- 
houses can be stated in the form of a standard mathematical 
programming problem as: 


Find the design vector X = \ . 




X 


(3.49) 


which minimizes the merit function f(X) 


subjected to the constraints g.(X) <_0, j = 1, 2,...,m 

J 


3.10.1 Objective function 

The nature of a thermal system design problem is such 
that there will usually be many designs that perform the 
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specified functional purposes adequately. The objective or 
merit function in a thermal design optimization problem repre- 
sents a basis for choice between alternate acceptable designs. 

The minimization of cooling and/or heating load and the minimi- 
zation of initial, plus operating and maintenance costs are 
taken as the objectives in the present work. In some cases, 
an objective function which is a linear combination of load 
(cooling and/or heating) and cost (initial plus operating and 
maintenance) may be considered for minimization. Whenever 
the cost is treated as an objective function, the economics 
model used in the estimation of costs plays an important role. 

In this work, an attempt is made to find the influence of the 
economics model on the optimum design of refrigerated warehouses. 

3.10.2 Design requirements 

The design restrictions are the upper and lower bounds 
on the insulation thickness and the bounds on the building 
envelope . 

3.10.3 Design variables- * 

The objective and the design requirements being known, 
the problem of optimum design of a refrigerated warehouse or an 
air-conditioned building can be cast as a mathematical progra- 
mming problem once the design variables are chosen. The 
thicknesses of walls, roof and insulation and the length, width 
and height of the warehouse are treated as design variables in 
this work. 
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3.11 Design Examples 

A computer programme has been developed for the 
automated optimum design of refrigerated warehouses by including 
all the factors stated in sections 3.1 to 3.8. The programme 
is quite general and can be used to optimize any refrigerated 
warehouse having any number of apartments . 

The optimum design of the refrigerated v^arehouse shown 
in Figure 3 .1(a) is considered to illustrate the effectiveness 
of the optimization procedure developed . The warehouse (with 8 
apartments inside) is assumed to have a known volume with all 
the four walls and roof exposed to atmosphere. The walls and 
roof are supposed to be constructed of bricks and reinforced 
concrete respectively with a 12 mm plaster on each side. All the 
walls and roof are insulated from inside. Two type of commonly 
used insulation, namely, expanded polystyrene and mineral wool 
are considered. The thermal properties of different materials 
and other design data are given in Table 3-2. The various 
aspects pertaining to the deterministic design of refrigerated 
warehouses are studied in six groups. The first four groups 
are devoted to parametric study. In the first group of problems, 
the thicknesses of roof and insulation are taken as design 
variables. Two types of insulation materials , as discussed 
earlier, are studied from economical point of view. The 
economical insulation material so found is used in subsequent 

design problems. In group two, thsree brick thicknesses are' 

' ■ 
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considered for the insulated walls. The strength of walls or 
roof is not considered | the design is based only on thermal 
considerations. In the third group, the absorptance of the 
exposed surface is considered as a design -variable although it 
could not be related to the cost in the absence of availability 
of cost data. After selecting the optimal thicknesses of walls 
and roof and the type of insulation, the warehouse orientation 
was changed (three different positions are chosen) for the same 
location in group four. Once the economical orientation is 
chosen, the optimal warehouse envelope is determined for a given 
volume in group five. The dimensions of the warehouse are taken 
from an. existing one at Kanpur, India. Finally in group six, 
a realistic problem is solved with all the internal loads . In 
all the above problems, the merit function is taken as the total 
cost, i.e., the sum of initial, running and maintenance costs. 
Another problem in group six is solved by considering the total 
cooling load as the objective for minimization. The statement 
and results of various optimization problems are given below . 

3.11.1 Group I 

The determination of optimum thicknesses of R.C.C. slab 
and its insulation is considered . The problem is solved for 
two different types of insulation. Two economics models, namely 
the annual cost with compound interest and the annual cost with 
simple interest, are used for each problem. Each of the 
problems can be stated as follows: 
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such that f(X) = f 1 (X) + f g (X) + f (X) + f 4 (X) 


subjected to ($) 
g 2 (X) 
g 3 ( x) 


= X t - 0.2 < 0 

= 0.10' - x 1 < 0 

= X 2 - 0.2 1 0 

= 0.005 - X 2 < 0 


minimum 


where X^ = R.C.G. slab thickness, X 2 = insulation thickness for 
roof (flat), f -j (X) = equipment cost, f 2 (X) = insulation cost, 

*4 - 4 - 

f^(X) = operating and maintenance cost, and f^ (X) = cost of 
R.C.G. slab. The costs of R.C.C., insulation and equipment are 
given in Table 3.2. These costs are converted to annual costs 
using Eqs. (3.43) and (3.46) for compound and simple interests 
respectively. The details of the four problems solved are given 
in Table 3.3(a) and the results of optimization in Table 3.3(b). 

It can be seen that in all the cases, the lower bound 
on concrete thickness is active at the optimum point. The two 
economics models give different results indicating their 
inherent characteristics. Since the optimum total cost with 
mineral wool as insulating material is lower than that with 
expanded polystyrene, mineral wool is recommended for insulation. 
The progress of optimization procedure is shown in Figure 3.3 
for problem 1 . 
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3.11 .2 Group II 

In this group three problems, each with a different 
brick thickness, are solved with mineral wool as the insulating 
material. Two different economics modes, namely the annual cost 
and the present worth, are considered for each problem. The 
optimization problems are formulated with insulation thickness 
as design variable as follows % 

Find X = {X^ 

such that f(X) = f 1 (X) + fg(X) + f^(X) + f ^ minimum 

- 4 * 

subjected to g^ (X) = X^ - 0.2 ± 0 

g 2 (X) = 0.005 - X 1 <_ 0 

where X^ = insulation thickness, f ^ (X) = equipment cost, 

-4 - 4 * 

f 2 (X) = insulation cost, f^(X) = operating and maintenance cost 
and f^_ = cost of bricks. The costs are given in Table 3-2. 

These costs are converted to annual cost and present -worth using 
Eqs. (3.43) and (3.42) respectively. The characteristics of the 
optimization problems are summarized in Table 3.4(a) and the 
results of optimization are shown in Table 3.4(b). 

Although the cost and energy requirements are minimum 
with 22.5 cm brick wall, due to the large envelope of the ware- 
house, the brick thickness of 33.75 cm is recommended for 
practical wall construction. This brick thickness corresponds 
to an increase of 4 % in annual cost and 5$ in energy requirement 



while the 45*0 cm brick thickness mdtes a 1 % increase in annual 
cost and hardly 1 % increase in energy requirement. Thus, if 
minimization of energy is taken as design criterion, then 
probably 45 .0 cm brick thickness might be more suitable than 
33.75 cm thick brick wall. The progress of optimization proce- 
dure for problem 9 is shown in Figure 3.4. 


3.11 -3 Group III 


A problem with insulation thickness and absorptance of 
exterior surface of the wall as design variables is solved in 
this group to assess the effect of absorptance of the surface. 
The brick thickness of the wall is taken as 33.75 cm and the 
insulation material is assumed as mineral wool. The annual cost 
is taken as the economics model. This problem is formulated 
as follows: 


Find X = 


4 * 

such that f(X) 


*- 





= f^X) + f 2 (X) + f 3 (X) 


minimum 


= X 1 - 0.2 -10 


• subjected to (X) 
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where = insulation thickness , X^ = value of absorptance , 

4 - 4 - 

f ^ (X) = equipment cost, f 2 (X) .= insulation cost, and f^(X) = 
operating and maintenance cost. These costs are converted to 
the annual cost by using Eq. (3.43). It is to be noted that the 
cost associated with absorptance (Xg) is not included in the 
above objective function. 

Since the lower bound on absorptance became active at 
the optimal point, it appears that lower values of absorptance 
would be more economical. However, the cost of achieving and 
maintaining any specified value of absorptance cannot be 
estimated with the available data. The progress of optimization 
procedure is shown in Figure 3.5. 

3.11.4 Group IY 

In this group the effect of orientation of the building 
is studied. The warehouse is assumed to be in three different 
orientations at the same location as shown in Figure 3.2(b) and 
the optimal cost corresponding to each orientation (to maintain 
the desired inside conditions) is obtained. Mineral wool is 
considered as the insulating material and the annual cost 
economics model is used. The optimization problem is formulated 
as : 
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such that f(X) = f i (X) + f (X) + f_(X) -»■ minim-urn 

I 2 5 

subjected to 0.005 f_X^ <_0.2 

0.005 <X 2 < 0.2 
0.005 < X 3 10,2 
0.005 1X 4 10.2 
0.005 1 X 5 10.2 

where X^ = insulation thickness for roof, X^ , . X = 
insulation thicknesses for the four walls, f ^ (X) = equipment 
cost, fglX) = insulation cost and f^(X) = operating plus mainte- 
nance cost. The result of optimization are shown in Table 3*6. 

By using the optimal roof and wall thicknesses obtained 
in earlier problems , the effect of orientation of the warehouse 
at a particular location is considered. The three orientations 
considered are shown in Figure 3.2(b). It is found from optimal 
design that the total cost is minimum for orientation 3, while 

the energy requirements are minimum for orientation Though 

% 

the cost reduction is hardly 0.6% for orientation 3 compared to 
orientation 2, the energy requirement is reduced by 1.6% for 

if ' ■ . 

orientation 2 compared to that of orientation 3. Thus orienta- 
tion 2 is preferred for the warehouse at the given location. 

The progress of optimization path for problem 12 is shown in 
Figure 3.6. 



3.11.5 Group V 


The determination of optimum envelope of the warehouse, 
i.e., finding optimal length, width and height of the warehouse, 
for a given volume and insulation thicknesses for roof and walls 
is considered in this group with the minimization of total cost 
as the objective. Only the cooling load due to building struc- 
ture is considered. The annual cost model is used for the 
investment policy and mineral wool is considered as insulation 
material. The mathematical formulation of the problem is as 
follows % 




X, 


X, 


Find X = J X 


X r 


X, 


X 


7 


such that f(X) =■ f ^ (X) + fg(X) + f^(X) 


minimum 


0.005 < X 1 < 0.2 


subjected to 
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0.6 < X g <_ 1 .2 

0.3 I x 7 1 0.6 

(0.455/Xg * X ? ) - 0.15 < 0 

0.1 - (0.455A 6 * X ? ) <0 

where X^ = insulation thickness for roof, X^, . Xg = insula- 
tion thicknesses for four walls , Xg = length of the warehouse , 

-K "K 

Zrj = width of the warehouse, f^ (X) = equipment cost, f^CX) = 
insulation cost and fg(X) = operating and maintenance cost. The 
design variables Xg and X^ are normalized and the volume of the 
warehouse is assumed to be 455,000 m . The results of optimi- 
zation are given in Table 3.7. 

It is found that, for a given volume of the warehouse, 
it is economical to have the maximum permissible height. The 
aspect ratio at the optimum point has been found to be appro- 
ximately 1.5. The energy requirements are reduced with the 
optimal envelope dimensions and insulation thicknesses by 10% 
while the aspect ratio changed from approximately 2.0 to 1.5 
and the height from 13.0 m to 14.9 m. The progress in the 
optimization procedure is shown in figure 3.7. 

3.11.6 Group VI 

In this group a realistic problem with all internal and 
external cooling loads is solved with two different merit 
functions. In the first case, the building envelope parameters 
and insulation thicknesses for walls and roof are optimized for 
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minimum total cost while in the second case, the same problem 
is solved for minimum total cooling load. The annual cost is 
taken as the economics model for the investment policy. The 
insulating material is assumed to be mineral wool. The mathe- 
matical formulation of the problem 17 is as follows: 

X 2 | 

. x 3 

Find X = X 

i 

X 5 
X 6 

X 7 

such that f(X) = L c minimum, subjected to 




g 2 5 

0.005 

^ x i 

< 

0. 

15 

§3 » 

S 4 : 

0.005 

i X 2 

< 

0. 

15 


§6 ; 

0.005 

v! 

< 

0. 

15 

Sy » 

§8 : 

0.005 

X 

V 1 

< 

0. 

15 


§io : 

0.005 

i X 5 

< 

0 . 

15 

Sir 

’§12 : 

0.6 < 

X 6 1 

1 . 

.2 


§13 

’§14 : 

0.3 1 

V I 

c— 

X 

0, 

.6 


§15 

♦ 

• 

(0.455/Xg * 

- x 7 ) 

- 0.15 <_ 0 

§16 

5 

0.1 - 

(0.455/X 6 

* X ? ) <0 


where L is the total cooling load and all other symbols have 
c 

the same meaning as in Group V, The results of optimization 
are given in Table 3.8. 
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In problem 16, the constraint corresponding to the 
upper bound on height became active while in problem 17, the 
upper bounds on the thickness of roof and two walls became 
active. The optimum point in problem 17 corresponds to an 
increase of 5.15% in cost while the cooling load is reduced by 
7.30%. In this case, the aspect ratio changed from 1.50 to 1.75 
while the height reduced from 14.9 m to 13.88 m. The progress 
in optimization procedure for the problem 16 is shown in Figure 
3.8. 

3.12 Kuhn-Tucker Conditions 

As stated in Chapter 2, for a constrained minimum, the 
Kuhn-Tucker conditions given by Eqs. (2.17) and (2.18) should be 
satisfied, let 

E = -vf(X ) 

[D] = [D ± . ] = (-jjJ) 

where i = 1 , . . . , n, 

3 e J, and 

J is the set of indices of active constraints 
then the value of ^ can be obtained as 

f = ([D] T [l])" 1 [l]' 1 E (3.50) 
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To illustrate that the final design points X obtained 
in different problems are optimum, Kuhn-Tucker conditions are 
applied to problem 17 (group YI) in which constraints 1, 5 
and 9 are active. For this problem S and [ D ] are found to be 

14.35 ^ 

7.82 

! 

12.80 j 
E _ } 8.44 V 

- T| 

33.46 
4.42 
3.70 | 

V 


i o o 

0 0 0 

0 1 0 

and [ D j = I 0 0 0 

0 0 1 

0 0 0 

0 0 0 [ 



The vector % given by Eq. (3.50) is computed as 


r 

F 

r 

1 


14.35 | 

j i 

■ ? 

> = ^ 

12.80 v 

S 

s 

33.46 

L j> 

r 
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Since all 5^ are positive, Kuhn-Tucker conditions are satisfied 
-*•* 

and hence X is guaranteed to he a relative minimum. 

3.13 Sensitivity Analysis 

In practice, a designer would he interested in knowing 
how the total cost or cooling load varies with a change in the 
design parameters. This type of sensitivity analysis will help 
the designer in manipulating the desin variables to suit some 
'i specific requirements. Further in some cases, the results 

obtained from the optimization procedure may have to he rounded- 
off to the nearest practical values of the design variables. 
Hence a sensitivity analysis of the cooling load/total cost with 
respect to different design variables is conducted at the 
optimum point. Thus the reference design is taken as the 
optimal design, and the values of the design parameters are 
changed by 50 % in steps of 10% on the positive and negative 
sides. The results of the sensitivity analysis for problem 17 
are shown in Figures 3.9 and 3.10. It can be seen from these 
figures that the cooling load is most sensitive to changes in 
the insulation thickness of the roof (X^ ) as well as the width 
of the building (X^) and least to changes in thickness of 
insulation of one of the walls (X^). The total cost is most 
sensitive with respect to the length of the warehouse (Xg) and 
least with respect to X^. 
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DETAILS OE PROBLEB, IS CONSIDERED IN GROUP I 
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OPTIMIZATION RESULTS OP THE PROBLEMS CONSIDERED IN GROUP VI 
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CHAPTER 4 


OPTIMUM DESIGN OP AIR-CONDITIONED BUILDINGS 

Air-conditioning involves the simultaneous control 
of temperature, humidity, air motion, air distribution, dust, 
bacteria, odour and toxic gases. Of these, the temperature, 
humidity and air motion are the most important ones. This 
implies that air-conditioning requires both cooling as well 
as heating. A building is a system which provides people 
with a comfortable environment. Human comfort is influenced 
by the psychological as well as physiological characteristics 
of human being-. There is no precise method of stating 
which thermal environmental conditions will affect the 
comfort feeling of a human being. The most important 
variables which influence thermal comfort are: (i) activity 
level (heat production in the body), (ii) thermal resistance 
of the clothing (clo-value), (iii) air temperature, (iv) 
mean radiant temperature, (v) relative air velocity, and 
(vi) water vapour pressure in ambient air. ASHRAE has 
studied human reactions to environment and gave the concept 
of an empirical index called "effective temperature" , which 
is an index expressing the combined effect of all the factors 
mentioned above on human comfort. 
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Fanger [57] has generalized the physiological basis 
of comfort so that for any activity, comfort can be predicted 
analytically in terms of the environmental parameters. 

Fanger 1 s comfort equation is of the form: 



(4.3) 
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where 


M = metabolic rate 

Aj^ = DuBois area (the surface area of the nude body 

(m 2 ) ) 


n = mechanical efficiency 
p = pressure of water vapour in ambient air 

Si 

t = air temperature (°C) 

Si 

f 1 = the ratio of the surface area of the clothed body 
to the surface area of the nude body 
t c ^ = mean temperature of outer surface of clothed body 
(°C) 


"^mrt = mean ra di«nt temperature (°C) 

O 

h = convective heat transfer coefficient (kcal/m hr °C) 
V = relative air velocity (m/sec) and 
I cl = thermal resistance of the clothing. 

Using the comfort Eq. ■ (4.1), it is possible for any type of 
clothing (clo) and type of activity to calculate all reasonable 
combinations of air temperature, air humidity, mean radiant 
temperature and relative velocity which will create optimal 
thermal comfort for persons in steady state conditions. The 
solution of Eq. (4.1) is also available in the form of comfort 
diagrams . 

However ASHRAE comfort standard 55-66 states thermal 
comfort as the condition of mind which expresses satisfaction 
with the thermal environment. It recommends the following 


conditions: 
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dry bulb temperature 
relative humidity 
air velocity 


= 73 - 77 °3? 

= 20 to 60 % 

.= 45 fpm maximum and 10 fpm minimum 


^mrt 


= 70 to 80 °E 


The common practice is to design comfort air- 
conditioning for summer to provide space conditions of about 
24 °C dry bulb temperature and 50% relative humidity. In fact 
Bevins and Gagge [58] have given new comfort charts which are ' 
primarily the combination of dry bulb temperature and water 
vapour pressure and Gagge and ITishi [ 59 ] have given physical 
indices for- thermal environment. 

The inside design conditions are dictated by comfort 
charts while the outside conditions are taken as stated earlier 
in Chapter 3. 


4*1 Computation of Cooling load 

• 

The various contributions to the cooling load are: 

(i) load from solar radiation, temperature difference of 
glass areas, exterior walls and roof. 

(ii) load due to heat gain through interior partitions, 

ceilings and floors of adjoining unconditioned space, 

(iii) load due to heat sources within the conditioned space 
such as occupants, lighting, power, equipment and 
appliances . 
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(iv) load from infiltration of outdoor air. 

(v) load from outdoor ventilation air. 

(vi) Miscellaneous heat loads. 

The method of computing heat gain from exterior walls 
and roof is same as the one given in Appendix A. The heat gain 
through the floor is neglected in air-conditioned buildings jj5(T[. 
The heat gain due to contributions (iii) , (iv) , (v) and (vi) is 
computed as described in Chapter 3. The heat gain due to glass 
areas is discussed in section 4.1.1. The total heat gain is 
converted into the cooling load by using the transfer function 
method . While evaluating the heating load , the light and human 
loads are not considered. Although theoretically the heating 
load would be reduced when the lights are switched on (with 
some time lag), this reduction is not considered in the present 
work. The plant is assumed to run during the period of 
occupancy only. 

4.1.1 Computation of solar heat gain through glass [1 ] 

The heat gain through glass due to the solar radiation 
incident on the outer surface and the temperature difference 
between outdoors and indoors constitutes upto 25$ of the total 
equipment load in the air-conditioning of many modern buildings. 
A procedure of estimating this component of load is presented 
in this section. Figure 4.1 shows the disposal of a quantity 
of monochromatic direct solar radiation 1^ incident upon a 
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single sheet of glass of thickness 1. A part of the incident 
radiation is reflected from the front surface and the remaining 
part is absorbed by the glass material. Because of successive 
internal reflections, the reflected, absorbed and transmitted 
radiation quantities are given by sums of infinite series. The 
total monochromatic transmissivity B^ is given as 



(1 - R f ) 2 a 

“ 12 2 

1 - R^a 


(4.4) 


where R^. denotes the fraction of each component reflected and 
a indicates the fraction of each component available for absor- 
ptions. Similarly the total monochromatic reflectivity is 
given as 


p x = 


R f + 


R f (1 - R f ) 2 a 2 

__ 

1 - 


(4.5) 


As <*.= 1 - B, - p. , Eqs. (4.4) and (4.5) can be used to obtain 


= 1 - R f - 


(1 - R f ) a 
1 - R f a 


(4.6) 


In order to evaluate the absorption coefficient, a, used in Eqs. 
( 4 . 4 ), ( 4 . 5 ) and ( 4 . 6 ), it is assumed that the absorbed radiation 
is proportional to the intensity of incident radiation and to 
the length of the path of the reflected beam. Thus 



K I dL' 
c x 


(4.7) 
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and 


a 



-K L 1 
. c 


(4.8) 


where K q is called the extinction coefficient.. 


Since 


cosO 


0 -£ i£^) V2 

n 


(4.9) 


where 0 is the angle of incidence of sun's rays, and n is the 
index of refraction of the glass, further the component refle- 
ctivity E^ may he found as [ 61 ] 




sin 2 (£ - &' ) 
sin 2 (6- + O' ) 


tan 2 ( 0 — O ' ) 
tan 2 (G + O’ ) 


J 


( 4 . 10 ) 


The use of two separated sheets of glass is common in buildings. 
Parmelee [62] has shown that for double glass, 

B B 

B = 1 - X1 — — (4.11) 

^’ 2 " P X1 P X2 


and 


P 

XI ,2 


P 

XI 


+ 



P 


X 2 


P 

XI X2 


( 4 . 12 ) 


where the subscript X 1 ,2 refers to the double glass combination, 
,X 1 refers to the first sheet of glass considered separately, 
and X2 refers to the second sheet of glass considered 
separately. 
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In heat transmission computations, it is more useful 

to know the solar absorption in each glass sheet rather than 

the absorption for both sheets. Thus we find 

[ 1 - (p + B )][1 - p ( p - B ) j 
XI XI X2 XI XI 


a. 


X,1 of 2 


1 - p p 

XI X2 


( 4 . 13 ) 


and 


°x,2 of 2 




P )] B 
X2 XI 


( 4 . 14 ) 


Once B and a are computed, the heat gain or loss Q by the 

o 

interior of the structure through the glass can be computed as 


Q = F B,I, + B^ + B-L + h. (t . - t. ) 

g sdd D D R R g,i i' 


( 4 . 15 ) 


Eq. ( 4 . 15 ) can be simplified further as 


(E a , I , + + ot—I— ) 


r o ■ 


+ 


(4.16) 


where 


U = 


1 


& + ( k ] 


(4.17) 


E = sunlit fraction of the window, dimensionless, 
s 

B^,B^ and B^ = transmissivity of direct, diffuse and reflected 
solar radiation respectively, 

Ig = incidence of reflected solar radiation 
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Ijj = incidence of diffuse sky radiation 

1^ = incidence of direct solar radiation upon a surface and 
a d J a D an< ^ °r. = a ^ sor P’ti v i^y of direct, diffuse and reflected 
radiations respectively. 

The heat gain Q is computed from Eq. (4-16). The data for 

O 

direct and diffuse solar radiations is used as available from 
the meteorological department and reflected radiation 1^ is 
assumed to be negligible. 

4.2 Reduction of Cooling Load 

The cooling load can be reduced by a proper design of 
the building by providing solar shadings, reducing air leakage 
and providing proper insulation to exterior walls and roof. 

The shadings can be designed either external or internal so 
that the solar radiations are reduced to a- minimum. The other 
internal loads like lights, fans and equipment can also be 
reduced by economizing their use and placing them at suitable 
places. However, for existing buildings, Guy [63] has 
suggested the following plan of action for optimizing the 
performance: 

(i) Reduce the primary loads on building systems 

(ii) Employ the most energy effective system predictable for 
the application 

(iii) Maintain the system for peak effectiveness 
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' (iv) Operate the system for minimum time periods 

(v) Educate tenants and operating personnel on energy 
conservation programme 

(vi) Retain the services of a professional engineer to survey 

the system and supervise the execution of a plan of action. 

4.2.1 Solar shadings 

The cooling load can he reduced hy designing the solar 
shadings. It can. he done in the following ways. 

i) Solar shading hy structural design 

The shading on glass windows can he reduced hy forming 
a canopy over the window head or hy placing vertical fins at 
the sides of the window. Some shading can also he achieved hy 
providing windows at suitable places in the v/alls. The surface 
finish, position and colour of the shading projection can affect 
its effectiveness. A high solar reflectance is desirable hut 
the position of the surfaces must be such that the radiation is 
not reflected hack on the glazing, 

ii) Solar shading hy blinds, louvres or drapes 

The direct solar radiation may also he excluded hy 
devices such as Venetian, canopy or roller blinds which can he 
positioned to cut-off the direct radiation. In addition the 
sky and ground diffuse radian can he reduced to some extent 
depending on the design of blinds. 
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iii) Shading by surrounding buildings 

The designer cannot control the type of surrounding 
buildings. However, if there are relatively tall buildings near 
enough to cast a shadow on to the building designed, there can 
be a considerable saving in the solar heat load. 

4.3 Hofee Control in Air-Conditioned Buildings 

The creation of an accoustical environment is as 
important as thermal environment for comfort in any air-condi- 
tioned building design. The designer must not only establish 
heating, cooling, and air flow performance requirements, but 
also accoustic performance limits for the various components 
and must specify pertinent installation details. It is also 
possible to isolate the building structure from vibrational 
forces by proper installation of equipment, ducts, etc. In 
fact, vibrations radiate in the form of noise when the frequency 
exceeds that of the lower limit of hearing (25 cps). It is 
possible to reduce the forces transmitted to the structure to 
a small fraction of the unbalanced force (of equipment) by 
using mounts having a sufficiently large static deflection. 

Indoor shading devices , particularly draperies', are 
effective in absorbing sound originated within the conditioned 
space. They have little or no effect in preventing outdoor 
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sounds from entering the room, lor excessive internally gener- 
ated sound, the usual remedy is to apply accoustical treatment 
to the ceiling and sometimes to add carpets as well. 

4 .4 Statement of the Problem 

In the present work the automated optimal design of 
an office building whose details are given in ASHRAE fundamentals 
[l5l, is considered for illustration. However, it is assumed 
that apart from the structural load, other loads are also brought 
to a minimum level. Two problems one with the objective fun- 
ction as the total cost of the plant (initial running and 
maintenance costs) plus the cost of thermal insulation, and the 
second with the objective function as total load are solved. 

The load is taken as the weighted average of the cooling ard 
heating loads. The weights are taken to be proportional to the 
number of months during which cooling or heating is required in 
a year. Whenever the average daily solar -air temperature falls 
below the desired inside temperature, then it is assumed that 
the building is to be heated in that month. The design day for 
such a month is chosen based on the minimum average solar-air 
temper ature and the hourly temperature and solar radian values 
corresponding to this date of the month are taken for computa- 
tions. The design day for the other months are chosen based 
on maximum value of the solar-air temperature as discussed in 
Chapter 3. The characteristics of the design days computed are 
given in Table 4.1 for Hew Delhi, India. 
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The plan of the office building considered for optimi- 
zation is given in Figure 4.2. The other details are given 
below: 

Occupancy : 85 office workers from 8 a.m. to 5 p.m. 

Lights : 17,500 watts, fluorescent and 4000 watts 

tungsten 

Four motors: 7.5 horse-power 

Inside conditions: temperature 24 °C and relative 
humidity 50$ 

Outside conditions: as given in Table 4.1 
Roof : flat, 0.1 m R.C.C. slab 

Floor construction: 0.1 m concrete on ground 
Fenestration: 0.85 m x 1 .5 m windows of regular plate 
glass with light colour Venetian blinds 
Front doors: Two, 0.75 m x 2.1 m (wood panels) 

Side doors : Two, 0.75 m x 2.1 m (wood panels) 

Rear doors : Two, 0.75 m x 2.1 m (wood panels). 

4.5 Problem Formulation and Numerical Results 

Two design problems, one with merit function as 
minimization of weighted combination of heating and cooling 
loads and the other with merit function as the total cost, are 
formulated and solved as nonlinear programming problems. 
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Problem 1 : Find the insulation thicknesses for the walls and 

roof of the office building for minimum load when upper and 

lower bounds on insulation thickness are prescribed. 

The load is computed as the weighted average of cooling 

and heating loads.' The weights are evaluated as follows: 

Let the solar-air temperature remain below the desired 

inside temperature for m^ months and above the desired inside 

temperature for the remaining (12 - m^ ) months in a year. 

Further let the heating load corresponding to each of the m^ 

months be q^_ and the cooling load corresponding to each of the 

(12 - m 1 ) months be q . Then the weighted average load 
i c ± 

considered (L ) is given by 


m 


i i. v + (12 - V I. v 


1 

r 


12 


(4-. 18) 


Mathematically the design problem can be stated as 

r 


Find it = < 


L 1 


X, 


X. 




such that f (X ) = L^ minimum 

subjected to g^gg: 0.005 5. 

g 3 ,g 4 : 0.005 iX 2 < 0.10 
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gpj, ggS 0.005 < X 3 <_ 0.10 
g ? , gg : 0.005 <X 4 £ 0.10 
gg, g 1Q J 0.005 <_Xg <_0.10 

where X^ is the thickness of roof insulation and X^, X^ , X^ and 

X c are the thicknesses of insulation of walls and 1 is the 
5 r 

weighted average of cooling and heating loads . 

The results of optimization are shown in Table 4.2. 

Problem 2 ; The second problem can be stated as follows ; 



such that f(X) = f 1 (X) + f 2 (X) + f^CX) minimum 

subjected to the same constraints as given in problem 1. 
where f^(X) = equipment cost, fg(X) = insulation cost and f^(X) 
= operating plus maintenance cost. The optimization results of 
this problem are given in Table 4 . 2 . 

It is interesting to note from Table 4 .2 that the 
insulation thicknesses of walls and roof at the optimum point 
are different in the two problems. In problem 1 (for minimi-* 
zation of load), maximum insulation thickness is obtained 


107 


for the roof indicating that cooling load has been dominant. 

In problem 2 (for minimization of cost), the roof has thinner 
insulation than the three exposed walls which indicates that 
the energy transfer during heating has been dominant. However 
the percentage cost and load reductions are less because the 
structural cooling and heating loads form only about 20 % of the 
total load. Depending on the design criterion, the solution 
of the first or second problem can be taken as the optimum 
solution. The progress in optimization procedure is shown in Figs. 
4-3 and 4.4. However, the results obtained in the present work 
are not compared with the ones given in reference [15] as the 
outside conditions are different in the present calculations. 

4.6 Kuhn-Tucker Conditions 


The Kuhn-Tucker conditions are applied to problem 1 
at the final design point X at which the constraints 1 and 7 
are active. The value of the vector I” is obtained by using Eq. 
(3.50). For this problem the vector E and matrix [D] are 
obtained as 


E 


2 


11 .83 
17.95 
13.63 
21 .72 
16.49 



t 
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r 


and [ D ] = 


1 

0 

0 

0 

0 


0 

0 

0 

1 

0 


J 


The vector 'f is computed as 


5 = 


f* 


r ^ 

{ 1 1 


— i 

00 

j 

= < 


! ^7 ! 

j 

21 .72 I 

L > 



> 0 


Since both the components ox f are positive, Kuhn-lucker 

J 

conditions are satisfied and hence the final design vector X 
is guaranteed to be a relative minimum. 


4.7 Sensitivity Analysis 

Figures 4.5 aid 4.6 show the effect of perturbing the 
design variables on the positive and negative sides about the 
opt imum point. It can be seen that the cooling load as well 
as the total cost is most sensitive to the insulation thickness 
corresponding to the roof (X^ ) and least to the insulation 
thickness of one of the walls (X^). 


DESIGN DAYS FOR AIR-CONDITIONED BUILDING DESIGN FOR NEW DELHI, 1967 * 
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CHAPTER 5 


FORMULATION AND SOLUTION OP PROBABILISTIC DESIGN PROBLEM 

Meteorological conditions play an important role in 
the design of heating or cooling systems. These conditions 
are Beyond the control of the designer and are I'andom in 
nature. A design which takes into account these random 
outside conditions and generate the desired inside conditions 
optimally has to be produced. Frequently some of the design 
parameters and properties of the materials used will also be 
random. Under these conditions, the deterministic approach 
for design will not be rational and a probabilistic design 
approach would be more suitable and realistic. Toussaint [ 64 J 
has evaluated the economic loss due to discomfort caused by 
random fluctuations in temperature, humidity, air velocity, 
noise etc. in air-conditioning systems. Green and Smith [21 ] 
have discussed the use of probability theory in the design 
of buildings. However, the probability of occurrence of a 
specified load value, which varies with a variation in the 
outside conditions, is not computed in all the above works. 

In this chapter the following problems are solved 
using the probabilistic design methodology : 

1 . Minimize the mean value of cooling load subjected to (a) 
upper and lower bounds on insulation thicknesses, and (b) 
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an upper bound on the standard deviation of the cooling 
load where the upper bound is taken as a fraction of the 
average value of the cooling load. 

2. Minimize the mean value plus three times the standard 
deviation of the cooling load subjected .to upper and lower 
bounds on insulation thicknesses. 

3. Minimize the mean value of cooling load subjected to (a) 
upper and lower bounds on the insulation thicknesses, (b) 

a constraint that the standard deviation of the load should 
be less than some constant times the average value of the 
cooling load and (c) a restriction that the probability 
of realizing the maximum inside temperature less than a 
given value should be more than a specified value. 

In these problems the probability characteristics 
of the cooling load which play an important role in probab- 
ilistic design have to be computed . The probability density 
and distribution functions of the cooling load can be found 
once the distribution functions of the meteorological and 
other parameters influencing the cooling load are known. 
Although this is possible theoretically, it poses a difficult 
problem in practice as the expression for the cooling load is 
quite complicated and the exact distributions of the random 
variables affecting the cooling load will not be known. 

Hence an approximate method, known as partial derivative 
method which makes use of only the mean values and standard 
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deviations of the independent random variables, is used in 
computing the probability characteristics of the cooling load 
in this work. 

5.1 De s ign Day 

The design day for the probabilistic design can be 
chosen based on a different criterion. It can be taken as the 
day of the month in which the value of mean plus three standard 
deviations of solar-air temperature has a maximum value. Twelve 
such days, each corresponding to a different month of the year, 
can be taken for design computations. 

A second definition of the design day can be taken as 
the fictitious day of the month which has the hourly temperature 
and solar radiation values same as the average hourly temper- 
ature and solar radiations of all the days of that particular 
month. The standard deviations of hourly temperature and 
solar radiation can also be computed for each ’day if needed. 

The first definition is used in selecting the design day in 
this work. The characteristics of design days computed are 
given in Table 5.1 for Dew Delhi, India. 

5.2 Probabilistic Description of Meteorological Data 

The meteorological data for a number of years is 
required to compute the mean and standard deviations of solar- 
air temperature . In the absence of availability of hourly 
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data for a number of years, the values available for one year 
are used for the mean values and standard deviations are 
assumed to be certain percentage of the mean values. 

Once the design days are chosen, the cooling load is 
to be computed on these days since the final load is taken as 
the weighted average of the individual loads corresponding to 
different months. The material properties and the building 
envelope parameters are also assumed as random variables. 


5 .3 Calculation of the Cooling load 


If all the parameters are deterministic, the heat 

gain can be calculated by the classical Threlkeld method D3 

according to the flow chart given in Figure A. 2 of Appendix A. 

The same procedure can be adopted for calculating the mean 

value of the heat gain also provided all the random parameters 

are taken equal to their respective mean values. However for 

computing the standard deviation of heat gain, the theory for 

finding the variance of a function of several random variables 

has to be used. According to the partial derivative method 

[22] used in this work, if F is a function of the n random 

variables X. , ..., X , then the standard deviation of F, a , 

1 tl Jj 

is given by 


9F 


F 


3X, 


>,2 2 r i F I 

— — _ y °v + v air J j w y 

y y Y -A- -1 o Y Y Y n 

*1 ’2 » * . • ’n - 2 1 X 1 ,x 2 ,... 2 


\2 _2 

- ) 


+ ... + ( 


3F 

3X 


n 


\2 2 

) 

X 1 ,X 2’.“’ X n n 


(5.1) 
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where 


__ is the partial derivative of the 


i 1 X 1 ,X 2 ,...,X n 

function F with respect to X^ evaluated at the mean values X_. 


of X. and a„ is the standard deviation of X. (i = 1, 2, 


» n) 


5.3*1 Mean and standard deviations of cooling load 


The expressions for the mean and standard deviations 

of cooling load can he derived as follows: 

The probability characteristics (mean and standard 

deviation) of solar-air temperature (t ) are given by 

0 



(5.2) 

(5.3) 


where a bar represents the mean value and a sigma represents the 
standard deviation of the parameter and the terms in brackets 
are evaluated at the mean values of the variables. 

The mean and standard deviations of solar radiation 
incident on the wall or roof per unit area I can be obtained by 
using the expressions given in section 3-3-3. The declination 
, total and diffused solar radiations (1^ and 1^) and absor- 
ptivity a are assumed to be random variables. The mean and 
standard deviations of the direction cosines (cos Z,cos'W and 
cos S) of the solar beam are 
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z = 

cos V sin 1 sin 6 + cos L cos 6 cos H ] 

(5.4) 

2 

CT Z " 

p 

(sin 1 cos 6 - cos L sin 6 cos H)^ 0 

(1 - cos 2 Z) 

(5.5) 

w = 

cos - ^ [ cos 5 sin H ] 

(5.6) 

2 

a W = 

sm 5 sm h _2 

9 Ts 

(1 - cos^ W) 

(5.7) 

s = 

cos - ^ [ 1 - cos 2 W - cos 2 Z ]^ 2 

(5.8) 

2 

a s = 

(cos 2 f sin 2 W + cos 2 Z sin 2 Z a 2 ) 

(1 - cos 2 S) cos 2 S 

(5.9) 

The mean and standard deviations of the incident angle 

9 can be 

expressed 

33 


9 = 

cos ^ [1 cos Z + m cos W + n cos 6 j 

(5.10) 


(l 2 sin 2 Z 


2 2 2 — 
+ m sin W 


W 


p 9 

+ n sin 5 




9 


(1 - cos 2 9) 


(5.11) 


Thus the mean and standard deviations of the incident radiation 
I are given by 

I = [(I T -T d ) cos 9 + I D ] I (5.12) 

a 2 = a 2 cos 2 9 a 2 + (l-cos9) 2 a 2 + [ (Tm-I^) cos& 

1 1 T D 


(5.13) 
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The mean, and standard deviation of the overall heat transfer 
coefficient, U, are given by 

1 


U = 


1 - + ^ + l 


fcu K 


h 


( E“ + I + iT } 4 [ ( V °h + (j 2 ) 0 i + ( ~4 } 

h ± k h Q h 4 n ± h K + - 


U 


•<k> 

O 


(5.14) 


( 5 . 15 ) 


For a composite wall, let suffix i represent the 
inside layer, hl, , m 2 , the first, second, n middle 

layer, and o the outside layer. Then the mean and standard 
deviations of (S) in terms of (S) of individual layers can 

jx 6 Q[ ^ 

be expressed as 

(I) + 0 + 0 + 0 + 


( 1 ) 


K eq 


In . , fl IN -...+(£) 

K m 


(5 - 16 ) 




eq 


K i 

K o 

K m^ 

K mg 

(\) 
n . 

*rH 

CNI Hi 
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,L 2 s 

+ (~t) 

K 4 . 

4 + (\) 

K i IT 
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(N) 


*<$ 
TT ‘ m 

_2 

0 K + *•' 


n 


g L ■ r V4 ; K 
o 0 K o 0 


(V 4 + 

K m n 


IN _ 2 
( N4 } K 
^ m n m n 


(5.17) 
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The mean and standard deviations of the coefficient t ^ in 

e ,m 

terms of solar -air temperature at different hours are given as 


t 


e ,m 


1_ 

24 




+ t 


e24 


(5.18) 


2 

a t 

e ,m 


1 _ 

24 



+ cr, 


"e2 


\ 3 

x e24 


(5.19) 


The mean and standard deviation of the product (KpC) g ^ and 
(S i) are given as 

££ fcj 


(K PC ) 


eq 


(h 

K 


3— [1.1(5) (K P C ) . + 1.1(5) (KpC; + ... + 

) K i K m„ 1 


eq 


(KpC) 


i.k5) (kpc)] + ^[(5) -o.k5) 

t r t t r t r 


K 


m 


h fls 


n 


(=) 

K 


K 


o 


K 


m. 


eq 


o.i(5) 3 

K 


(5.20) 


m 


n 


(KPC) 


1 


eq 


_ . {1.1 (IpC) + 1.1 (LpC) + ... + 1.1 (LpC) + 

1 n 

* .q 


(Kp°) „[(=) - 4)=c0-1- ••• - 0.1 (i) ]} a K + 

vJ T7* 'Tr t r l \ 


K o K m 


K 


1 


hl 


n 




eq 


— 5 0-1 PC) 

(1.1 PCT a 2 + [ m 1 

(I, 2 11 & 

K eq 


K 


eq 
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(KPC) „ , 2 , 

— — 2 (£*!)] <2 + ... +[ 

(5) V m 1 

K 1 


(1.1 PC) 


m 


n 


eq 

O.I(KPC) 


(5) 

K 


eq 


„ t ,2 9 • ( PC ) 2 ? 

— 2 tC) J 4 +...+[ — — 2 J o| + 

(i) V m n (5) 


K 


n 


eq 


K 


(KPC) 2 V 2 


0 1 rrs A — +...+[ 


[ — — ] (0.1 v. 


eq 

(KPC ) 2 


(5) 

K 


K -mu 

1 


eq 


(5) 

K 


] 


eq 


m n 2 o • ^ • "Ik. ) p o —2 2 \ 

(0.1 z- 2 ) <y| 4 _ - (c? a ; + P? <£ ) H 

K V, ,3^ 2 

m„ n (— ) 


1 p i 1 C i 


n 


K 


eq 


(1 * 1L m } 2 

r Hi /~' 2 - 2 2 rt 2 ^ 

[— J ( c m a + P™ % ) + 


(5) 

K 


'm w C 

n m n n m n 


eq 


(Kj 


°' {(z) - O.l(z) - ... - 0. 1 (z) }] 


fhs K K 

(z) 0 m 


K 


K 


1 


m 


n 


eq 


(Sf % + 7 o °C ) +[^4-^<(=) - 0.1(5) 

° P ° ° °° (5) K o K m 


1 


o.k5) > 

K »n 


(PC) T 2 o 

(=> ] 4 

(5) K o 
K e4 


( 5 . 21 ) 
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CO— 1/2 T 

( S l>e 9 “ & 


( 5 - 22 ) 


eq 


2 

0 (S L) 
v n 'eq 


s (KPC) a Z T . +i(l) » (£p0)-' 


(-) 8 I 

V eq 

eq 


J n^ PW eq °(KPC) 


eq 


(5.23) 


The mean, and standard deviations of the coefficients S^K are 
given as 


(S„K)„„ = (t 3 (K-P0) eq ) 1/2 


n 'eq 


(5.24) 


S n K 


= ? a (K PC) 


(5.25) 


eq 


The mean and standard deviations of the coefficients Y^, 

and can be evaluated in a similar way . 

The mean and standard deviations of the coefficients 

M , 1ST - and t are given as 
n n s jXi 

M n = ^ ( \l cos “l T + *** + *e24 C0S “l x) + (T e1 C0S W 2 T 


+ . . . + t g2 ^ cos OSgx) 


(5.26) 


M 


n 


1 2 2 

(cos to- T + COS C0 o t) CT, + . • • + (cos 03 - T + 

12 1 2 X e1 


\2 2 
cos 03 2 t) cr t 


(5.27) 


e24 


U. 


n “ 12 


(t^^ sin co^t ■+■ •••“!" t^ 2 ^ ® in o 3 .j ,x) (bg-j sin co^t 


+ . . . + t e24 _ sin « 2 t) 


(5.28) 
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n 


1 2 2 
j 2 (sin y + sin t) a ^ 


el 


+ ... + ( s in a) ^ t + 


2 2 

sin u t ) c.|. 


e24 


e ,n 


(M 2 + N 2 ) 

n n 


1/2 


(5.29) 


(5.30) 


2 

a t 


e ,n 


— _ -1 _ 

(M 2 + N 2 ) (M 2 . al 

v n n ^ n ] 


+ N 


■ m 2 


n 


n * N 


aS ) 

n 


(5.31 ) 


The mean and standard deviations of the angular displacement 


0 , angle X^ and the decrement factor 4 n are given by 


2 


n 


-1 z 

tan ' 1 -S 


n 


(5.32) 


a 0. 


Y 2 

n 


n 


Y 2 Y 2 + z 2 ' 
n n n 


y2 O n-2 

rt 2 , / n / In a 2 

z ° ^ ^ 


(5.33) 


n 


= tan 


-i N 

-1 n 


M 


n 


(5.34) 


^n 


M 2 0 N 2 

( S ) a*~ + / n — \ a 

M 2 + l 2 N n M 2 + l 2 ~"n 

n n n n 


2 

M 


(5.35) 


n 


n 

2 


U 


"n 


V 2 

n CT 2 + V CT 2 

U 2 V n U 4 U 


(5.36) 


(5.37) 
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Finally the mean and standard deviation of the heat transfer 


q^ are given as 


l'i = 11 ‘C t e,m + l 1 Vi cos( “l T -^1 * V + 


X 2 \ 2 C0S ^“9 T - $9 - *p) + • • • ] - ^ } ( 5 * 38 ) 


*i 


lt e,m + X 1 Vi 00s( “l T - F 1 - V + x 2 V 2 oos( “2 T 


-V z - * 2 ) + ... - t ± > «2 + D 2 a 2 


U 


e ,m 


2 t 2 1 cos 2 (o) 1 t - ~ ♦ 1 ) a 2 x + U 2 X 2 cos 2 (w 1 t - 

^ - V a t . +,?2 *1 T e,i sin2( "i T ~ F i • V 

e » ! 


r 2 . 2 - 

[ °0 + 0 *, 3 + 


( 5 . 39 ) 


Since the heat transfer rate through the total surface area 


Ajjis 


Q. = A.q. 
1 11 


the mean and standard deviations of Q. are given by 


( 5 . 40 ) 


Q ± = Vk 


( 5 . 41 ) 


i ■ \ + ^ i 


(5.42) 
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To evaluate the mean and standard deviations of t m , 
M n and N , the integral is first approximated by summation over 
all - the discrete intervals of time and then the partial 
derivative rule is applied. However, these quantities can 
also be evaluated directly as follows: 

The expression for M is given as 


M 


n 


12 


24 

} t g cos to n T dT 


By letting X = with the upper limit of integration as y , 


the mean value X can be exnressed as 


X = M = Y 2 / t g cos co^t dT = t e ^ cos “n T dr 
x o o 




sin a y 


(5.43) 


^2 


The expected value of (X - X) can be written as 


E [(X - X) 2 ] = E[{— ^ t e C0S t °n T dT “ 12 ^ t e cos “n T dx} 3 

o o 



{t ( T ) - t } { t ( 6 ) - t } ] (5.44) 

Hence the standard deviation of X, can be obtained as 

of = 777 / Y di / Y dB cos “n T cos “n e K t t (5.45) 

A 144 o o e’ e 

where k, . ( T , 8 ) is covariance function of t g at different 

times t and & . * * 



131 


If t (x) and t (g) are independent raniom variables, Eq. (5-45) 
becomes 

2 rY ff o 

a X = 144 ' dT / cos 0> n T cos u n e a t S(t - g) (5.46) 

where 6 is the Dirac . delta which is equal to 1 if x = g 

and 0 if x ^ g . Hence 

2 "| 2 f* y 

• ’ • a v = TTT / dT cos u„x ex. / \ j dg cos M a fi.( t - g) 

X 144 o n V T ' o n 

(5.47) 

Thus cr^. can be obtained by numerical integration and the 

continuous values of t needed can be interpolated from the 

© 

hourly values available. However, it requires a large computer 
time to evaluate the integral a number of times. Hence the 
integral involved is approximated by a summation (over hourly 
intervals) for calculating the mean and standard deviation of X. 

5.3.2 Internal loads 

The mean values of internal loads are evaluated as 
discussed in Chapter 3 and their standard deviation is assumed 
as 5% of the mean value. However, if the data is available, 
then the mean and standard deviations can be computed. 

5.4 Problem Formulation and Numerical Results 

In order to illustrate the applicability and effect- 
iveness of the probability based design procedure in the 
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optimum design of thermal systems, three problems are formu- 
lated and solved. All these problems correspond to the design 
of refrigerated warehouses whose specifications are given in 

Chapter 3* 

Problem 1 : Find the insulation thicknesses for the walls and 
roof for the warehouse for minimum cooling load, when the 
upper and lower bounds on the insulation thickness are pres- 
cribed and the standard deviation of the load is constrained 
to be less than a fraction (k) of its mean value. Mathema- 
tically the problem can be stated as a standard optimization 

problem as': 

->• 

Find X = 

such that f(X) = L - minimum 

subjected to 0.005 £^i i 0.2 

0.005 1 x 2 - °* 2 
0.005 < X 3 <0.2 
0.005 1 X 4 1 0.2 
0.005 1 X 5 i 0.2 

°Tj ' “c i ° 




X. 


X, 


X, 


X, 


x r 


J 
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where X^ = insulation thickness for the roof, 

X^jX^jX^jX^ = insulation thicknesses for the walls, 

L c = mean value of the cooling load, 

= standard deviation of the cooling load, and 
c 

k = a constant . 

The results of optimization obtained by using mineral wool as 
the insulating material are shown in Table 5.2. 

The same problem is solved with the upper bound on 
the insulation- thicknesses changed to 0.15 m (instead of 0.20 
m) , and the results of optimization are shown in Table 5.3* 

It can be seen that the cooling load is reduced by 
18.27% in the first case and 15-90% in the second case. The 
total annual cost in the first case is 4.6% higher than that 
of the second case while the load requirement is less by 3.3% 
in the first case. Hence, depending on the situation, either 
of the two solutions can be chosen. The progress of optimi- 
zation procedure is shown in Figure 5.1. 

Problem 2 : Find the insulation thicknesses such that the mean 
value of the load plus three times its standard deviation is 
minimized, when upper and lower bounds on the insulation 
thicknesses are prescribed. The problem can be stated as a 
standard optimization problem as 
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Find X = 


4 x 


3 1 


h 

X(- 
V_ 5. 


such that f(X) 


(l c + 2>n: 0 ) 


minimum 


subjected to 0.005 £X^ £0.2 

0.005 £ X 2 <_ 0.2 
0.005 < X ^ < 0.2 
0.005 < X 4 < 0.2 
0.005 < X 5 < 0.2 


where X^ = insulation thickness for the roof, 

X 2 ,X^,X 4 ,X^ = insulation thicknesses for the walls, 

L = mean value of the cooling load, 

W 

a T = standard deviation of the cooling load, and 
^c 

k = constant, taken as 0.02 in the present case. 

Here also mineral wool is taken as insulation material. The 
results of optimization are shown in Table 5 .4 and Figure 5.2. 

It is obvious from this figure that maximum drop in 
the objective function value as well as in the penalty function 
value takes place in first iteration. In subsequent iterations, 
this drop, particularly in the case of objective function, is 
relatively small. It was observed that the insulation thickness 
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gradually increases. But the effect of this increase, after a 
particular combination, is not very marked in reducing the load; 
however this increases the cost more. 

Problem 3 s Find the insulation thicknesses such that the mean 
value of the cooling load plus three times its standard 
deviation is minimized, when upper and lower bounds are pres- 
cribed on the insulation thicknesses and the probability of 
temperature of the inside surface of the wall less than a 
prescribed value is required to be greater than a specified 
value. This problem can be stated as follows: 



such that f(X) = (L c ) + minimum 
subjected to g^ GO = - 0.15 i 0 


g 2 CO 

= 0.005 - X 1 1 0 

g 3 (x) 

= X 2 - 0.15 lo 

s 4 Cx) 

= 0.005 - X 2 < 0 

g 5 (X) 

= ■ X, - 0.15 1 o 

5 
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g 6 m 

= 0.005 - X < 0 

5 

g 7 (x) 

= X 4 - 0.15 10 

s 8 ($) 

= 0.005 - X^ < 0 

e g m 

= X c - 0.15 < 0 

3 ” 


= -0.005 - X < 0 
5 

s-n(x) 

= P(tg <_t) ip 

g 12 0O 

= 0 T - kL 1 0 

P. c 


where 

and X^ = insulation thicknesses corresponding 
to roof and four walls respectively, 

1 = mean value of the cooling load, 

o 

0 , = standard deviation of the cooling load, 

^0 

k = a constant , 

tg = inside maximum surface temperature of any wall 
or roof 

t = prescribed temperature value, and 

p = prescribed value of probability. 

The probability constraint can be converted into 
an equivalent deterministic form as suggested by Hati and Rao 
[65]. Let the probability constraint be written in general 
form as 

s ■ ■ - 

P [ g^'(x) < b. J > p. (5.48) 
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By defining a new random variable y . as 

J 


^3 = *3°° = - g 3 (X) 


(5-49) 


Eq. (5-48) can be written as 


J p(y 1 )dy i p i 

O J J 


(5.50) 


where p(y.) is the density function of the random variable y. 

3 . ^ 

By expanding E . (X ) about the mean values of X and retaining 
only the linear terms, one gets 


UP . 


y. = E.(X) - P.(X) + J sq ( j. , 


(X, - x ± ) 


(5.51 ) 


Here it is assumed that y . is a normally distributed random 

3 

variable and its mean and standard deviations are given by 


- *3® 


SP . 

^ ? ^x’i 
1 1 


2 2 -, 1 / 2 
X i 


(5.52) 

(5.53) 


with the transformation of variable 0 ^ and noting 

that 




_L_ „-tV2 


f e u ' - dt = 1 , 


(5.54) 


Eq. (5-50) reduces to 


1 2 * 

— e d0 > 




1 


-t 2 /2 


T}. (p ) 

a 3 


dt 


(5.55) 
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where ^ (p.) depends on the probability level p.. Thus 
3 D 


< l b (p . ) 

a y “ a V 

*3 ■ 


i.e. - y i - (p. ) o < 0 (5.56) 

j J 

Eq. (5.56) can be written as 

_ 3F. 2 0 1/2 

S-j( X ) - WE (I °i. ] s 5 (5.57) 

C X 1 A 1 


Thus the constraint reduces to the form of Eq. 

(5.57) where the value of $ (p.) is obtained from standard 

a 3 

probability (normal distribution) tables corresponding to the 
assigned probability level p . . 

For the inside surface temperature tg, first the wall 
and roof temperatures for same insulation thickness are 
computed and the one (walls or roof) having the maximum temper- 
ature at any time during 24 hours is picked up. Then the 
partial derivatives of this temperature with respect to each 
of the random variables, namely, the outside temperature, solar 
radiation incident and diffused, wind velocity (outside film 
coefficient) and inside film coefficient, are obtained numer- 
ically and the constraint value is evaluated using Eq. (5.57). 
This type of constraint can help in checking the instantaneous 
maximum inside surface temperature and hence the radiations 
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between wall and the stored commodity etc . at any instant which 
may affect the stored product adversely . 

By using mineral wool as insulating material, the 
optimization problem is solved. The results are given in 
Table 5.5 and the progress of optimization procedure is shown 
in Figure 5.3. The maximum inside surface temperature (t) was 
specified as 6.0°C and the probability of actual temperature 
exceeding 6.0°C was required to be not more than .3333%. The 
inside temperature is to be maintained as 0°C . Hence it is 
expected that there will not be any bad effect on the commodity 
stored in the absence of instantaneous cooling. The final 
maximum surface temperature (at the optimum point) came out to 
be 3.8°G. In fact this type of restriction may be quite useful 
for commodities which are very sensitive to temperature 
fluctuations. 

5.5 Kuhn-Tucker Conditions 

The Kuhn-Tucker conditions are applied to the final 
design point X of problem 3 in which the constraints 1 , 3 and 
7 are active. For this problem % and [d] are found to be 

^ 52.30 
13.70 

E = J 6.50 

12.23 
7.30 

. v 
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and [ D ] = 


1 

0 

0 

0 

0 


0 

1 

0 

0 

0 


0 

0 

0 

1 

0 




J 


The value of the vector % is computed, losing Eq. (3*50), as 


C = 


*1 

^3 f 
C 7 


52.30 
= i 13.70 


> 


-*■ 

0 


12.23 


V. 


Since all are positive, Kuhn-Tucker conditions axe satisfied 

-*■ * 

and hence the final design vector X is guaranteed to he a 
relative minimum. 


5.6 Sensitivity Analysis 

i 

Figures 5.4 and 5.5 show the effect of varying the 
design variables about the optimal point on the positive and 
negative sides by 50% in steps of 10%. It is observed that the 
cost as well as cooling load is most sensitive to the insula- 
tion thickness of the roof (X.) and least to the insulation 
thickness of one of the walls (X^ ) . 
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OPTIMIZATION RESULTS FOR PROBLEM 
WITH UPPER BOUND ON INSULATION THICKNESS 
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CHAPTER 6 


APPLICATION OP EXTREMAL DISTRIBUTIONS 

In most of the engineering design problems, design 
variables are random in nature. In many situations, concern 
often lies with the largest or smallest of a number of random 
variables . Success or failure of . a system may rest solely on 
its ability to function under the maximum demand (load) or 
minimum capacity (strength), not simply the typical values. 
Floods, winds, temperatures, solar radiations and floor 
loadings are all variables whose largest value in a sequence 
may be critical to several engineering systems. Similarly 
the smallest values of strength of materials will also be 
critical in some situations. A systematic study of the maximum 
or minimum values of the random variables influencing the 
design may lead to a better design. Further the concept of 
finite life design can also be incorporated and thus a better 
use of materials can be made and technological obsolescenece 
can be cared for. A basic argument against the study of these 
extreme values is that as their probabilities are small, the 
error involved in their computation may be veiy large . However 
it has been found by Gumbel [24] that extreme values are much 
more reliable than the median whose value increases with . * 
simple size. Although extremal distributions have been applied 



152 


to several engineering disciplines in the past, no application 
has been reported in the field of thermal ^rstems design. The 
application of extremal distributions in the design of heating 
and cooling systems is considered in this chapter. 

In the present work, the hourly data is used to fit 
three types of extremal (asymptotic) distributions for temper- 
ature and solar radiation. Since type III distribution for the 
largest -values has been found to fit the daily and yearly 
maximum values of dry bulb temperature and daily maximum values 
of solar radiation very well, it has been used in the optimi- 
zation work. The optimum design of a refrigerated warehouse 
is considered with the mean value of the cooling load as the 
objective and an upper bound on the probability of occurrence 
of a specified extreme (permissible) value of temperature inside 
the warehouse. Upper and lower bounds on the design variables 
as well as an upper bound on the standard deviation of the 
cooling load are also considered in the formulation of this 
problems . 

6.1 Extremal Distributions 

Let X be a random variable and n tests be performed 
to find the probability characteristics of X. Let the results 
of the test (values realized for X) be arranged in increasing 
order of magnitude X.j , X£, ..., so that the smallest value 
realized for X be X.j and the largest X^. Let these n tests be 
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repeated several times and the results X JJ , X!| , X” ' , . .., 

XI, X", X 0 , . . . , X ' , X" , X" 1 , ... he fitted into a distribution 
F x (x). Once the distribution of X is known, the distributions 
of the smallest value X^ and the largest value X^ can also be 
obtained. The essential conditions are: (i) Statistical 
variates are dealt with (ii) The parent distribution [3? x (x)] 
from which the extreme values have been drawn and its para- 
meters remain constant from one sample to next (iii) The 
observed extremes should be extremes of sanples of independent 
data. 

Distribution of the smallest value X^: 

P(X^ > x) = P(all n of the X^ > x) 

If X^ are independent , 

P(X^ > x) = P(X^ > x) P(X 2 > x) ... P(X^ >x ) 

= [1 ~ E x (x)][1 - (x) ] ...[1-P x (x)j 

In the specified case where all the X^ are identically distri- 
buted with cumulative distribution function F x (x), 

P(X 1 > x ) = [1 - F x (x) j n (6.1) 

F x ^(x) = P(X 1 < x) = 1 - [1 -F x (x)] n (6.2) 

Distribution of largest value X^: 

As in the case of the smallest value X^ , if X^ are 
assumed to be independent and identically distributed, 
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P(X n <x) = P x (x) = [F x (x)f (6.3) 

Knowing from past experience the distribution of X^ (say, 
temperature) in any (i^* 1 ) year, one might need the distribution 
of, say, XcjQ, the largest temperature in 50 years, this being 
the design life time of a proposed thermal system. Equations 
(6.2) and (6.3) represent the exact probability distributions of 
the extreme -values in terms of the parent distribution function, 
E x (x ) . 

If the conditions of independence and common distri- 
bution hold among the X^, then in most of the practical cases 
the shape of the distribution of extreme values (X^ or X ) is 
relatively insensitive to the exact shape of the distribution 
of X^. In such cases limiting forms (as n grows) of the distri- 
butions of X^ and X^ can be found, which can be expected to 
describe the behaviour of that random variable even when the 
exact shape of the distribution of X is not known precisely. 
These limiting distributions are called asymptotic distributions 
The theory of asymptotic distributions may also be applied to 
many practical situations even when the assumption of independ- 
ence is not strictly satisfied provided sufficiently large 
samples are used for obtaining the parameters of the extremal 
(asymptotic) distributions. 

following are the specific distributions which are 
applied widely in engineering applications. 
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6.1.1 Type I distribution for largest value 

If the parent distribution is of the type 

F x (x) = [ 1 - e“ s ^ ] (6.4) 


where g(x) is an increasing function of x, then the distribution 
function of the largest value, Iy(y), is given by 


I'y(y) = I’y (y) = exp {- exp [- a 1 (y-u)*3 } , -» <_ y < « 

(6.5) 

and the probability density function by 

- a. (y-u) 

fy(y) = a -| ex P [- a-, (y-u) - e ] (6.6) 


The parameters and u in Eqs. (6.5) and (6.6) are to be 
estimated from the observed data. 


6.1.2 Type I distribution for smallest value 

% 

The probability distribution and density functions for 
the smallest value are given as follows 

1*2 (z) = ]? x (z) = 1 - exp {-exp [- a^(z-u)]> , -« 1 z 1 00 

1 ' (6.7) 

and 

f 2 (z ) = exp {^(z-u) - exp [ o 1 (z-u)]> (6.8) 

6.1.3 Type II distribution for largest value 

If the parent distribution is 
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1 

F x (x) = 1 - g(^) x > 0 . (6.9) 

then the probability distribution and density functions of the 
largest value (Y) are given as' 

K 1 

ly(y) = \ (y) = exp [- (“) ] ( 6 . 10 ) 

n ^ 

and 

K K^+1 

f Y (y) = U (S) exp[- (|) ] y>0 (6.11) 

where the parameters and u are to be determined from the 
observed data. 


6.1.4 Type II distribution for smallest value 

The probability distribution and density functions for 
the smallest value are given as follows 


P Z (z) 

and 

f z ( Z ) 


1 - exp [- (|) J 

TT K 1+ 1 



( 6 . 12 ) 

(6.13) 


6.1.5 


F x (l) 


Type III distribution for largest value 

If the parent distribution is of the form 
K i 

= 1 - c(o> 1 - x) , x i ; K 1 > 0 


(6.14) 


then the probability distribution and density functions are 
given as 


a ) 1 - y K i 

\(y) = exp E~( — - -) ] , y i w 1 


(6.15) 


and 
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fy (y ) = 


K 


1 


0) 


1 


“1 - U • 1 - U 


■-ir 


K 

oil - y i 

exp {-( — ) } , 

Ol-J u 


y i g>i 


(6.16) 


where the parameters and u are to be determined from observed 
data. 


6.1.6 Type III distribution for smallest value 


If 


* x (x) 


= c(x - e) 


¥L J 


(6.17) 


then 


K„ 


P z (z) 


V z) 


= 1 - exp t- ( " }, z <. 


u - e 


(6.18) 


and 


f z (z) 


K 


U - e 


1 £ 


V 

S'* , fZ - e 

-) exp {-(JJ-T7 


K 1 

) } , Z < e 


(6.19) 


where the parameters u, e and are determined from the 
observed data. 

These three types of distributions for largest and 
smallest values are used in different displines of engineering. 
Type I distributions are applied to describe strength of brittle 
materials and hydrological phenomena, type II distributions are 
used to represent annual maximum wind and several meteorological 
and hydrological phenomena while type III distribution is used 
to study strength of metals in tension and fatigue. 
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6.2 Pitting Extremal Distributions for Meteorological Data 

The above-mentioned distributions are fitted into 
maximum daily temperature and solar radiations for a year for 
two cities of India, They are also tried with yearly maximum 
temperature. The unknown parameters of the distributions are 
determined from the observed values. Por fitting the data, 
the expressions (6.5), (6.10) and (6.15) are rewritten as 
follows. 

The distribution function for type I (for largest value) 
is given by Eq. (6.5) as 

Py(y) = exp {- exp [- (y-u) ]> (6.20) 

Taking logarithm of Eq. (6.20), we obtain 

- In Py(y) = exp [- (y-u) J (5.21) 

i.e. In (jT-jyj) = ex P C- a 1 (y~ u )3 (6.22) 

By taking logarithms once again, we get 


—In ln( 
°1 



- y + u 


(6.23) 


i.e. y = - n m + u (6.24) 

Eq. (6.24) represents a straight line if y is plotted as 
ordinate against -In ln( p (y'J ^ 35 hbscissa, provided the parent 
distribution given by Eq. ( 6 . 4 ) holds true for the phenomenon. 
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Similarly the distribution function for type II (for 
largest value) is given by Eq. (6.10) as 



(6.25) 


By taking logarithms of Eq. (6.25) twice and rearranging the 
terms, we obtain 


In Y = 


— jr In ln(j, 


TyJ 


) + In u 


(6.26) 


1 

Thus if In y is plotted as ordinate against -In ln( |- ^y ) as 
abscissa, then Eq. (6.26) can be represented as a straight line 
provided the pa-rent distribution given by Eq. (6.10) holds for 
the ranlom variable. 

Lastly the distribution function for type III (for 
largest value) as represented by Eq. (6.15) is 


%(y) 


= exp 


, *1 ' l 

^ - u 



(6.27) 


where y <_ and > 0 


Again taking logarithm of Eq. (6.27) and rearranging the terms, 
the following expression is obtained . 


ln(o> 1 - y) = jj- In ln( ^-|y y) + ln(w ¥ 


u) 


(6.28) 


1 

Eq. (6.28) with ln( w-p y) as ordinate and In ln( p --- ^. -y) 85 
abscissa can be represented as a straight line if the parent 
distribution given by Eq. (6.15) holds true. 
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Prom the available hourly temperature and global solar 
radiation data, first the maximum of daily temperature and 
global solar radiation are selected. The maximum values so 

obtained are used to determine the values of abscissa and 

# 

ordinate at each point for the three types of distributions. 

For this the N individual observations (x^ , x^, . . . ,x m , . . . ,x^.) 
are arranged in increasing (decreasing) order of magnitude for 
maxima (minima). Then the N fractions y , m = 1, 2, . . . , N, 

are calculated and the points [{x m , -r^j -j}]are plotted on 
extremal probability paper. Once the different points are 
plotted, then a straight line is fitted through these points 
using least squares method. The error can also be computed 
between the exact point and the corresponding point on the 
straight line . The magnitude of this error is a gauge to 
measure which of the distributions fits best in the observed 
data. In Eq. (6.28), was assumed to be little higher than 
the observed maximum-most value. The results obtained are 
given in Tables (6.1) and (6,2) and are shown graphically In 
Figures 6.1 to 6.6 for the dry bulb temperature and global solar 
radiation of New Delhi and Roorkee, India. These results 
indicate that type III distribution fits best for temperature 
as well as solar radiation data. 

To find the difference between extremal distribution 
and the commonly used normal distribution, the daily maximum 
values of temperature are fitted into normal distribution and 
the resulting density function is compared with the type III 
extremal density function given by Eq. (6.16) in Figure 6.7. 
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6.3 Application to Design of Thermal Systems 

To demonstrate the use of extremal distributions in 
the design of thermal systems, the optimum design of a refri- 
gerated warehouse is considered. For this the type III extremal 
distribution based on yearly maximum values of dry bulb temper- 
ature for 25 years (for Poona, India) is used. All the three 
types of distribution are fitted to this data using Eqs. (6.24), 
(6.26) and (6.28), and the results are shown in Table 6.3 and 
Figures 6.8 to 6.11. It is found that type III distribution 
for the largest value fits the data best. Then the hourly 
values of temperature are computed using this maximum value. 

The temperature rarge and the percentage daily range are taken 
as the average of the values of twelve days of the year, one 
selected from each month on which the temperature is maximum 
(in that month). 

The hourly temperature is calculated as 

t, = 1; - DLY * PDLY h (6.29) 

h max n 

where t^ = temperature at any hour h, 

DLY = daily range of temperature , and 

PDLYj i = percentage daily range at h"^ hour of the day. 

The hourly values of solar radiations are taken over the day 
having the maximunvalue of mean plus three times the standard 
deviation. In this design procedure only one such critical 
day is taken as the design day. 
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6.4 Problem Formulation and Numerical Results 


Here the aim is to find insulation thicknesses ^ 

^ for 

the walls and roof and maximum design temperature for 

minimum mean cooling load, (ii) minimum total cost and (ij_i) 

minimum of the weighted sum of total co st and cooling l OQ ^ 

Further, upper and lower bounds are prescribed on the -s 

•insu- 
lation thicknesses, the probability of occurrence of the 

maximum design temperature is restricted to be less than a 

prescribed value and the standard deviation of the load i 

constrained to be less than a certain fraction of its me^n 

value. The problem with the first objective can be stated as 

a standard optimization problem as: 



such that F(X) = L minimum, subjected to 

v 


g v g 2 : 0.005 IX 1 < 0.15 

g 3 , g 4 : 0.005 iX 2 i 0.15 

g 5 , g g : 0.005 1 x 3 1 0.15 

gy» gg *• 0.005 1 £ 0.15 

g 9’ g l 0 : 0.005 <_ X 5 <,0.15 
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© 


where X^ is the thickness of roof insulation, Xg, X^, X^_ and 

X c are insulation thicknesses of walls and X c is the maximum 
5 6 

temperature, By (xg ) is the probability of occurrence of 

temperature less than Xg , is the maximum possible value of 

x 6 , Oy is the standard deviation of the cooling load, L c is 
c 

the mean value of the cooling load and K^, are constants. 

The value of the probability ly(xg) is calculated 
using Eq . (6.15) with the values of constants evaluated from 
the equation of straight line (fitted by least squares method) 
as follows. 

The equation of the straight line obtained is 

ln(w 1 - y) = 0.2992861 In In( gr^ y) + 1.231853 (6.31 ) 

By matching the coefficients of Eq. (6.31) with those of Eq. 
(6.28), we obtain 


A- = 0.2992861 or K. = 3 .3401 (6.32) 

K 1 1 

and InU-j-u) = 1 .231853 or (» 1 - u) = 3-4212 (6.33) 


© Xg represents the particular value taken by Xg. 
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The value of the constant in the constraint - Fy(xg) < 0 
is taken as 0.05. The three optimization problems stated 
above are solved using the procedure described in Chapter 2 
and the results are given in Table 6.4. The progress of 
optimization for the three problems is shown in Figures 6.11 
to 6.13. 

The minimum cooling load obtained in case of the 
first problem is 3,19,613.7 kcal/hr while it was 3,03,854.7 
kcal/hr in the case of probabilistic design (Chapter 5). Thus 
there is an increase of 5*08% in the cooling load due to a 
change in the design criterion. The three problems solved to 
illustrate the application of extremal distributions give 
different cooling loads and annual costs. The reduction in 
cooling load is 18.6% -in problem 1 while it is 12.7% and 
17*5% in the second and third problems respectively. Similarly, 
the cost reductions in the three problems are 2.7, 5.2 and 
3.6% respectively. It is interesting to note that in the 
first problem, the percentage load reduction is maximum while 
that of cost is minimum compared to the other two problems. 

In problem three, where equal weightages are given to cooling 
load and cost, the amount of reduction obtained in t be cooling 
load or cost lies in between the values obtained in the other 
two cases. This shows that an objective of the type used 
in problem three has to be used whenever both energy conser- 
vation and investment are equally important. 
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6.5 Alternate Formulation of Design Problems 

The problems of design of thermal systems can also be 
formulated in alternate forms using the concept of return 
period in extremal distributions. The return period T(x) is 
defined as: 

T(x) = 1 _ p( x ) 

X 

where ]^(x) is the distribution function of X. Physically the 
return period indicates the number of observations during which, 
on the average, there is one observation equalling or exceeding 
x. Its value increases with that of the variate (x). It can 
be seen that to eveiy distribution function there corresponds 
a return period and vice-versa. If the observations are made 
at constant intervals of time, the return period being a number 
of observations, is also a time, measured in the same units. 

The return period corresponding to any extreme value 
can be obtained from the probability paper as shown in Figure 
S.loox corresponding to any return period the extreme value of 
the parameter can be obtained. Now the design problem can be 
formulated with this extreme value and the related return 
period as the life of the system. As an example , if the return 
period corresponding to the extreme temperature is to be 
obtained, a horizontal line FE is drawn from the desired value 
of maximum temperature, F (43«6°C) and from the point of 
intersection E of lines FE and AB (line fitted by the least 
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square error method for type III distribution), a vertical line 
is drawn until it intersects the abscissa at the line G. The 
value of the return period corresponding to this point is 
found to be 20 years in Figure 6.8. Since this temperature 
(43 .6°C ) is expected to be exceeded only after twenty years, 
the life of the system can be taken as this return period , 

namely , 20 years* 


6.6 Kuhn-Tucker Conditions 

The K uhn -Tucker conditions are applied at the final 
design point X* of problem 3 in which constraints 1 and 14 are 
active. For this problem the vector E and matrix [D 3 are found 
to be 


and 



[p ] = 


r 


< 


V- 


17.51 
13.73 
11 .82 
9.62 
15.46 
18.75 


1 0 ^ 

0 0 

0 0 

0 0 

0 0 

i 

0 1 
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The value of the vector K is computed as 



17.51 

18.75 


> 5 


Since both are positive, Kuhn-Tucker conditions are satisfied 
and hence the final design vector X is guaranteed to be a 
relative minimum. 


6 .7 Sensitivity Analysis 

Figures 6.14 and 6.15 show the influence of changing 
the design variables (from the optimum point) on the total cost 
and cooling load. It is seen that the cooling load as well as 
the total cost is most sensitive to the outside maximum design 
temperature. Among the insulation thicknesses, the roof 
thickness (X^ ) has the maximum effect on the cooling load and 
total cost . 
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CHAPTER 7 


CONCLUSION AND RECOMMENDATIONS 


7 . 1 Conclusions 

(i) The optimum design of refrigerated warehouses and air- 

conditioned buildings using multidimensional optimization 
techniques has been considered in this work. The 
solution procedures to be adopted according to determin- 
istic and probabilistic design philosophy are indicated 
and illustrated . The designs are based on twelve 
critical days (instead of one used in conventional 
procedures ) selected on the basis ofi maximum average 
solar-air temperature . The minimization of cooling load 
as well as total cost are considered in obtaining the 
numerical results. In the computer programme developed, 
a linear combination of these two objectives is also 
considered so that proper weights can be assigned to the 
two objectives depending on the requirements. 

(ii) When the design of roof, which is a major source of heat 
transfer, is considered, it was observed that minimum 
possible thickness of roof with appropriate insulation 
thickness would be most economical. 

(iii) From the results of the problems in which the brick 
thickness is considered as a design variable, it has 
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been found that the optimum point corresponds to minimum 
thickness of the wall with proper insulation thicknesses. 
However, from physical (strength) considerations one may 
have to use larger wall thicknesses. 

(iv) The orientation of the building at any given location is 
not found to have an influence greater than 0.6% in cost 
and 1 .6% in cooling load at the optimum design. 

(v) When building envelope parameters are also in the control 
of the designer, it will be more economical to use an 
aspect ratio of 1.5 to 1.75 (1.5 from cost consideration 
and 1 .75 from energy consideration) and larger heights 
to reduce the cost as well as cooling load. It is 
interesting to note that this conclusion agrees with the 
present trend of the design of warehouses where larger 
heights are used. 

(vi) The effect of economics models on the minimum cost 
design has also been studied. Since different cost 
models influence the optimum design, one has to decide 
about the economics model before designing a refrigera- 
ted warehouse/air -conditioned building. 

(vii) The use of random parameters like geometric parameters, 
material properties and environmental conditions such as 
temperature, solar radiation and wind velocity, in the 
optimum design of refrigerated warehouses is illustrated . 
This procedure assumes that the cooling load follows 
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normal distribution, which is justified from the central 
limit theorem. Further, this assumption simplifies the 
calculations. The procedure developed, however, is quite 
general and can be used even when the cooling load 
follows a distribution other than normal distribution. 
Although the results obtained according to probabilistic 
procedure are not substantially different from the ones 
obtained by using deterministic procedure, the probab- 
ilistic procedure is expected to be more realistic and 
rational. 

% 

(viii) Three types of extremal (asymptotic) distributions are 
fitted into the maximum hourly temperature and solar 
radiation data and type III distribution (for largest 
values) has been found to be most suitable for describing 
these phenomena. An optimum design procedure using 
extremal distributions is also presented and illustrated 
by considering the design of a refrigerated warehouse by 
using the yearly maximum temperature data for 25 years. 

As expected, the optimum designs obtained by using 
extremal distribution have been found to be more conser- 
vative from the point of view of both cooling load and 
total cost compared to deterministic design and probab- 
ilistic design based on normal distribution. 

(ix) The use of a weighted combination of cooling and heating 
loads in the design of air-conditioned buildings is 
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demonstrated. It is found that the insulation thickness 
corresponding to roof is maximum when the total load is 
taken as the criterion function while it has a lesser 
value when cost is taken as the objective. 

(x) Whenever the stored commodity is more' sensitive to fluc- 
tuations in the inside conditions (and large penalty or 
risk is involved), the inside conditions can be restric- 
ted to be at the desired levels in probabilistic design 
procedures . 

(xi) The interior penalty function method has been found to 
be quite satisfactory for the optimum design of thermal 
systems. It has been observed that the largest decrease 
in the objective occurs in the very first iteration in 
all the cases. Hence if one is interested in finding 
only near-optimal solutions (rather than the exact 
optimum solutions), the process can be terminated after 
five or six iterations only thereby saving a substantial 
amount of computer time. 

(xii) The results of sensitivity analyses conducted about the 
optimum design points are expected to be useful whenever 
the designer is required to round-off the design vari- 
ables to the nearest available values. Further this 
information would be useful in identifying the lesser 
sensitive design variables so that the designer need not 
consider them while dealing with the design of similar 
systems . 
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7 .2 Recommendations for Future Work 

(i) The design of equipment* can also he considered along 
with that of building in the optimization procedure. 

(ii) The analysis procedure used in the calculation of cool- 
ing load can be refined by using a more accurate method 
like matrix method. 

(iii) The efficiency of other optimization methods like the 

method of feasible directions, geometric pro gr amma rig arid 
dynamic programming can be studied. 

(iv) The availability of discrete values of insulation thick- 
nesses can be directly considered by using integer non- 
linear programming techniques . 

(v) The optimum design of thermal systems using the concept 
of return period can also be studied. 

(vi) The concept of overall reliability of the building 

system can also be incorporated with the optimum design 
procedure. 

* 

(vii) The accoustical considerations can also be incorporated 
along with thermal considerations in the optimum design 
of air-conditioned buildings. 


Compressor, condenser, evaporator, ducts, piping etc. 
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APBBNDIX A 

HEAT TRANSFER THROUGH A HOMOGENEOUS WAIi 


In this appendix the periodic transfer of heat through a 
wall made up of a single homogeneous material is considered. 

Eigure A.1 shows the schematic problem. In this analysis it is 
assumed that (i) the heat flow is one dimensional (along x- 
direction) , (ii) the wall is homogeneous with constant material 
properties, (iii) the surface coefficients Ir and h Q are constant, 
(iv) the solar absorptivity of the outside surface is independent 
of the angle of incidence and is constant, (v) the variation of 
outside air temperature t Q and incident solar radiation I are 
periodic, and (vi) the internal thermal environment is constant. 

At any location within the wall, the one dimensional 
conduction equation can be written as 


3t 3 2 t 

~ Ct o 

3x^ 


(A.1) 


where 

t = wall temperature at time t and distance x, 

K 

a = thermal diffusivity = -j^, 

K = thermal conductivity of the wall material, 
p = wall density and 

G = specific heat for the wall material. 
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Equation (A.1) can be solved by using the two boundary 
conditions: (i) At the inside surface, 


*i = 


- K(— ) 
v ax 7 


X=Ij 


= h. (t . - t. ) 
i w,i i 7 


(A. 2) 


where h^ 


combined convection and radiation heat transfer 


coefficient for the inside surface of wall, t w ^ = temperature 
of the inside surface of thi 
and at the outside surface, 


of the inside surface of the wall and t^ = inside air temperature. 


% = -*<|E> X=0 - h o (t e - (A ' 3) 

where h = combined convection and radiation heat transfer 

o 

coefficient for the outside surface of the wall, t,„ = tempera- 

1 W jO 

ture of the outside surface of the wall and t = solar -air 
temperature given by 


t = t + (A. 4) 

e o h Q 

which can also be expressed as 


t = r« + ^ i cosCrn. t - ) + t 0 cos (« 0 t - 1> ) + ... 

e e,m e,i i l e ,2 2 2 

(A. 5) 

where t = 24-hour mean value of solar-air temperature and 
e,m 

t = harmonic coefficients (n= 1, 2, ...). The complete 
0 f n 

solution of Eq. (A.l) has been given by Alford, Ryan and Urban 
[3], The temperature of the inside -wall surface t w i may be 
expressed as 



20 i 


= t ± + £-[u(t 


~ lU(t e,m - t i ) + Ve ,1 cos( V ~ *1 “ + 


Y 2 t e 2 cos ( “2 T - ^2 “ $2) + ***^ 


(A. 6 ) 


where the overall heat transfer coefficient U is given by 


U = 


1_ + Ii + i- 

h. K h 
1 o 


(A.7) 


*e ,m “24 { t e d ' 


(A. 8) 


Since the values of t g can only be computed at hourly intervals 
from the available data, Eq. (A. 8) is taken as 


= 24 CVl + V + ••• + *e24 ] 


(A. 9) 


The expressions for V^, and are given by 


Vi 

V < Y n + Z n> 


(A. 10) 


where n = order of the harmonic in the Eourier series of Eq. 


(A. 6 ) , 


0i„ 1/2 

s n - 


(A. 11) 


where w n is the angular velocity of the sinusoidal wave. 
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Y n = 


h h. h h. 

(~^2^2 + O cos ^ ^ sin k S L + (— 0 d - 0 - 1) sin S L cos S I 
?^T n n 2S 2 K 2 n n 


2sjr 


and 


(h o + h i } 

q 77* cos S Ii cosh S Xj 
o ii xi n 

n 11 


(A. 12) 


n 


k h * h h . 

(~ "o 9 + 1 ) sin S I cosh Si - ( — 6 \ - 1 ) eos S L 
OC*^ 31 II o c* TT n 


XI' 


n 


h + h. 

sinh SJD + ( ° g v - ) sin Si sinh SI 


n 


S n K 


n 


n 


(A.13) 


0 n = tan 


-i z 

-1 n 


(A. 14 


n 


In Eq. (A. 14), sin$ n has the sign of and cos0 n has the sign 


Again 


where 


e ,n 




p.1/2 


o 9 »/ ^ 

= (BT + ]» ) 
n n 

(A. 15) 

1 24 

T2 S *6 C0S “n T dT 

0 

(A. 16) 

24 

U/t e Sin „ n T dT 

(A. 1?) 


= w/12 radians per hour or 15° per hour, and “ n = n “ • 

In Eqs. (A. 16) and (A. 17), the integrations can be 
replaced by s ummat ions because only hourly values can be 
computed for the solar-air temperature. Hence 
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M n ~ 12 ^ t e cos °n T 

0 

1 24 

S n = 12 I *e sln »n T 

1 N 

ip = tan - srj^ 
t el M 


(A. 18) 

(A. 19) 

(A. 20) 


In Eq. (A. 20), the quadrant in which ♦ lies is determined by 

the requirements that sin $ must have the sign of N and cos ^ 

n n n 

the sign of M^. 

The rate of heat transfer to the interior can be 
expressed as 


= VWV • (A - 21) 

By substituting Eq. (A. 6) for t y/ Eq. (A. 21) can be written as 


4 i = tJ{ t t e,m + Ve.l 00s( V “ *1 ‘ + X 2 t e,2 C “2 T ‘ *2 ~ 


where 


0 £ ) + * * * ] “ tj_ ) 

V 

_ n 
~ U 


(A. 22) 
(A. 23) 


Hepe is called the decrement factor and 0^ is the angular 
displacement or lag between a harmonic of the solar-air temper- 
ature aid the sane harmonic of the inside surface temperature. 
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Equivalent temperature difference: 

The form of Eq. (A. 22) is interesting. Although q^ 
may be continuously changing, it may be calculated by multi— 
plying the overall heat transfer coefficient U for steady state 
heat transmission by an equivalent temperature difference which 
accounts for the periodic variation of solar-air temperature t 
and the heat storage characteristics of the wall. Equation 
(A. 22) can be written as 

q^ = U A TgQ (A. 24) 

where 

4T EQ = [ %,m + Ve,1 oos( "l T ' +1 - + 

cos (u 2 t — 4*2 “ ^ 2 ^ + • • • J - t^ ( A. 25 ) 

where A Tj,q is known as the equivalent temperature difference. 

In the present wodc, Eq. (A. 22) is used to compute the 
heat gain through the homogeneous wall. The order of harmonics 
is taken as 2 because it was found that with 4 and 6 order 
harmonics, there is negligible difference in the heat transfer. 

The flow chart for the heat gain through exterior 
walls and roof along with the other subroutines used in the 
computer programme is given in Figure A. 2. 







r \ 

START y 


ead NCODE* , A ~L , V ,K , p ,C ,DATE , MONTH ,WALL 
ORIENTATION, (I , J) ,I D (I , J) ,h ± 

(I = 1, 24, J = 1, 1 2 )** 


Read t m£DC ,DLY,PDLY(I) , 
I = 1 , 24 


Read t Q (I,J) 
1 , 24; J = 


Compute t 0 (I,J-) = 
t max -DLY*PDLY(I) 


CALL SRDATA (Mouth, Date , Orientation) 
This subroutine returns the required' 
solar radiation I 


Compute ho using Eqs. (3.14) and (3.17) 


Compute t e using Eq. (A. 4)] 


Compute t. _ using Eq. (A. 9) 
6 9 111 
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*If NCODE = 1 Hourly outdoor air temperature data is available 

£ 1 Hourly temperature data is to be computed using 
maximum outdoor air temperature tniar ■> daily 
range DLY and percentage daily range EDLY 

**I = 1,...,24 (number of hours per day) 

J = 1,...,12 (number of months per year) 
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DETAILS OP OTHER SUBROUTINES USE IN THE COMPUTER PROGRAMME 


FTN 


PTM 

GLHT 

GRADN 


: To compute cooling load through different 

external and internal sources and objective and 
0-function 

: To evaluate maximum inside surface temperature 

: To evaluate the heat gain through glass 

: To evaluate the gradient of 0-function using 
forward difference formula 


MINIM4 : To implement the David on-Pletcher-Powell 
variable metric method of unconstrained 
minimization 


SL0PE2 

SRDATA 


STEPK3 


To evaluate the slope of 0-function 

To compute solar radiation incident on any 
wall or roof at any hour of the day 

To evaluate step length using cubic inter- 
polation technique of one dimensional 
minimization. 


FIGURE A. 2 PLOW CHART POR THE HEAT TRANSFER THROUGH 
WALLS AND ROOF 
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APPEMDIX B 

DATA USED IN COMPUTATIONS OP HEAT GAIN* 


TABLE B;i 

VALUES OP EQUATION OP TIME POE DIPPEBENT MONTHS 


Month/Date 


|Equation of time 
{ (min : sec . ) 


Month/Date 



! Equation of time 
{ (min : see . ) 



January 

21 

— 

11 


18 

Pebruary 

21 

— 

13 


28 

March 

21 

— 

7 


19 

April 

21 

+ 

0 


08 

May 

21 

+ 

3 


32 

June 

21 

- 

1 


48 


July 

21 

— 

6 


25 

August 

21 

- 

1 


18 


21 

+ 

7 


30 

0 ctober 

21 


15 


06 

November 

21 

+ 

13 


55 


21 

+ 

1 


32 


TABLE B.2 

VALUES OP DECLINATION POE DIPPEBENT MONTHS 


Month/Date \ 
\ 

L 


January 

21 

- 

20.00 

July 

21 

+ 

20.60 

Pebruary 

21 

- 

10.80 

August 

21 

+ 

12.30 

March 

21 


0.00 

September 

21 


0.00 

April 

21 

+ 

11 .60 

October 

21 

— 

10.50 

May 

21 

+ 

20.00 

November 

21 

- 

19.80 

June 

21 

+ 

23.45 

December 

21 

- 

23.45 


Month/Date 


Declination 
(degree) 


Declination 

(degree) 


*A11 the data given in this appendix is in P.P.S. units. Thi s 
data has been converted with M.K.S. units in actual computations. 
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TABLE B.3 

AVERAGE AIR CHANGES PER 24 HOURS FOR STORAGE ROOMS 
DUE TO DOOR OPENINGS AND INFILTRATION 


f 

Volume | 

cu . f t . | 

- T 

Number 
per 24 

of air changes 
hours 

200 


44.0 

300 


34.5 

400 


29.5 

500 


26.0 

600 


23.0 

800 


20.0 

1000 


17.5 

1500 


14-0 

2000 


12.0 

3000 


9.5 

4000 


8.2 

5000 


7.2 

6000 


6,5 

8000 


5.5 

10000 


4.9 

15000 


3-9 

20000 


3.5 

25000 


3.0 

30000 


2.7 

40000 


2.3 

50000 


2.0 

75000 


1 .6 

100000 


1.4 


Note: For heavy usage multiply the above values by 2.0. 
For long storage multiply the above values by 0.6 
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DATA OF POOD PRODUCT 
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APPENDIX C 

DATA BELATED TO TEE CONVERSION OP HEAT GAIN TO COOLING LOAD 


TABLE G.1 


CONVECTIVE AND RADIANT HEAT GAIN TO COOLING LOAD 


t 

Heat gain component }% 

- _ _ . f 

1 

radiant |% convective 

I 

}% latent 

t 

Heat transfer through trans- 
parent surfaces (without inside 
shading) by solar radiation 

100 

0 

0 

Heat transfer through trans- 
parent surfaces (with inside 
shading solar radiation) 

58 

42 

0 

Heat transfer through trans- 
parent surfaces (with or Without 
inside shading) by air to air 
convection and conduction 

0 

100 

0 

Heat addition by fluorescent 
lights 

50 

50 

0 

Heat addition by incandescent 
lights 

80 

20 

0 

Energy addition by people 

40 

20 

40 

Energy addition by machinery 
or appliances* 

Energy addition by infiltration 
and ventilation 

0 

Convective and latent 
account for 100% 

Heat transfer through interior 
partitions or exterior walls 
and roofs 

60 

40 

0 


*Percentage will vary according to the device being considered. 
The sensible portion could vary from 20 to 80% * radiant (and 
80 to 20% convective) depending on its surface temperature. 
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COEFFICIENTS OP ROOM TRANSFER FUNCTION 3 
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